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Linear Elastic Stability 


A Cnitical Analysis of Methods 
By Hans Ziecter, E.T.H., Zurich 
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IV. Buckling by Compression 


Straight rods and shafts present a variety of applications of the preceding — 
theory. 


19. Euler’s Buckling Cases 


A straight rod, subjected to a sufficiently small axial thrust P (Figure 19.1), 
has a stable equilibrium configuration in which it is slightly compressed. If its 


* 


Pp 
Z 
ARLES 
p 
x 
Figure 19,1 Figure 19.2 
Rod under axial thrust. Buckling cases. 


equilibrium is disturbed, it executes axial, torsional and flexural vibrations. As 
P increases, the equilibrium, for a certain value P,, becomes unstable, and the 
rod, being heavily bent by its load, buckles. 

The so-called buckling load P, depends on the constraints (Figure 19.2). 
The reactions, provided the constraints are rigid and frictionless, do no work. 
If damping, in accordance with Theorem 16c, is disregarded, the internal forces 
are conservative. So is the load, provided it is not only of a constant magnitude 
but also retains its axial direction independent of the shape of the deflection 
curve (Example 2.3). Hence, the problem is purely nongyroscopic. From 
Theorems 14 it follows that any load P => P, is critical and that P, can be 
obtained by either of the static methods. ' 

Let the rod be homogeneous and prismatic. Then, provided it is supported 
symmetrically with respect to the principal planes (x, z), (y, z) (Figure 19:39 
the displacements of the deflection curve in either of these planes are indepen- 
dent. Thus, we are confronted with the buckling problem of L. EULER [25]}), the 
solution of which can be obtained by individual treatment of either projection, 
while compression and torsion may be disregarded. The result is the well known 
formula of L. EULER 

Rahat, (19.1) 
where / denotes the length of the rod, f its smaller flexural rigidity and ka 
numerical factor taken from Table 19. 


1) Numbers in square brackets refer to the Bibliography, page 184. 
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Table 19 
Euler's buckling loads 
| Case | 1 2 3 | + | 5 | 


| Factor k | 4 | il | 0-25 | 2.045 | 1 | 


20. The Influence of Constraints 


By Theorems 14 (Table 14) any purely nongyroscopic buckling problem is 
characterized by a buckling load P, which represents the limit between the 
stable and the unstable regions and is obtained as the smallest value of P for 
which V is not positive definite. In order to distinguish it from other kinds of 
critical loads (particularly from the so-called buckling loads of higher order), 
we might call it a genuine buckling load. 

Let A and B be two purely nongyroscopic systems differing only by their 
constraints. If B is subjected to all the constraints of A and, in addition, toa 
number of workless constraints, A and B have the same potential energy. Any | 
configuration of B compatible with the constraints also is an admissible con- 
figuration of A (while the reverse is not true). Thus, if V in the system A is 
positive definite, so it is in the system B; hence, the buckling load of B can not 
be smaller than the one of A. We have, therefore, 

Theorem 20a. If, in a purely nongyroscopic system, workless constraints are 
added, the buckling load does not decrease (but generally increases). 

From Theorem 6 it follows that Theorem 20a, formulated for the smallest 
critical load, remains valid in the presence of gyroscopic and dissipative forces. 
Thus, we obtain 

Theorem 20b. If in a stationary, noncirculatory system workless constraints 
ave added, the smallest critical load does not decrease (but generally increases). 

The constraints illustrated by Figure 19.2 are workless, provided that friction 


is negligible, Their number increases from Case 3 to 1 and from Cases 3 and 5 
respectively to Case 4. Table 19 shows the corresponding increase of k. 


21. Asymmetric Constraints 


If the constraints (Figure 19.2) are asymmetric with respect to the prin- 
cipal planes of the rod, the projections of the deflection curve, being connected 
by the end conditions, can not be treated separately. There are various types 
of asymmetric constraints; perhaps the most important one is represented by 
a large bearing, the axis of which is rotated in the end section through an 


arbitrary angle # with respect to the y-axis. 
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In Figure 21.1 Case 5 is shown, both axes being rotated through the same 
angle # which can be restricted to the interval 0 < # < a/2. The buckling load [27} 
is given by (19.1) again, the numerical factor k (by which P, exceeds the EULER 
load of Case 5) depending on the angle #. In Figure 21.2 the square root of this 
factor is plotted versus #, the parameter A? = «/f denoting the ratio of the flexural 
rigidities. For #= 0, k =1; the rod buckles in the plane perpendicular to the axes 
of the bearings, F, being the Euler load in Case 5. If # increases from 0 to 2/2, the 
buckling load is raised by the increasing stiffness with respect to the axes of the 


19 Z 7 ve 
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, l!/// oh 
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13 /| 
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Figure 21.1 Figure 21.2 
Rod with rotated hinges. Buckling load for rotated hinges. 


bearings, the deflection curve generally being a space curve. The highest value 
k = 4, attained by rods of slender cross section (42 > 4) for = x/2, corresponds 
to buckling in the plane of the axes, P, being the Euler load in Case 1. 


Other examples of asymmetric constraints have been examined by P. Fit- 
LUNGER [28]. Like the case just treated, they are obtained from cases of 
Figure 19.2 by adding or removing workless constraints; hence, their buckling 


loads, according to theorem 20a, are limited by the corresponding values of 
Table 19. 


22. Twisted Rods 


If the rod is twisted in the unstrained state by a constant angle w per unit 


length, its displacements in either direction are connected by the differential 
equations. 


In Case 5, the buckling load is found [29] to be given by (19.1), the numerical 
factor k (by which P, exceeds the Euler load) depending on the total angle of 
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twist, ® = /. In Figure 22.1, k is plotted versus % for several values of A? = «/B. 
On the whole, the buckling load, owing to the stiffening effect of the greater 
flexural rigidity «, increases with the angle of twist, k tending, for # > co, to 
2 22/(1 + A*) = 2. For a given value of 42, the factor k is determined by an 
infinity of successive arcs. The reason is obvious: near # = 2 7 the buckling load 
of the first order (of the untwisted rod) begins to exceed the second order load, 
and since this process is repeated with higher orders, P, is successively determined 
by buckling loads of increasing order. i 


‘ 1 Fees 


vy | =e 
QO 45 90 135 100 205 270 35 30 YO5 450 495 HO 585 6H 675 720 
Figure 22.1 
Buckling load of a twisted rod. 


Case 3 has been solved by E. LUscuer [30]. Other cases of twisted rods 
have not been treated so far. 


23. Tangential Compression 


The rods considered in 19 to 22 are purely nongyroscopic. If, however, the 
thrust P, although of constant magnitude, has no constant direction but 
remains tangential to the deflection curve (Figure 23.1), the 

problem (Example 2.4) may be circulatory. 


I 
7 VA Resolving P as shown in Figure 23.1, we obtain, in the case 
of a prismatic rod with symmetric constraints, 
a DG 5 MES RE oie (23.1) 


While V is conservative, H represents the circulatory constituent 
of P. From Figure 19.2 it follows that H, apart from Case 3, 


ee | does no work. Hence, Cases 1, 2, 4, and 5 are characterized by 
Rod loaded genuine buckling loads F,, given by (19.1) in connection with 


tangentially. | Table 19. In Case 3, however, the existence of a genuine buckling 
load is uncertain. Besides, the critical loads according to Theo- 
rem 17, must be calculated by means of the kinetic method. Eventually, latent 
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instabilities must be taken into account, appearing, by Theorem 11a, whenever 
a nontrivial equilibrium position exists. 

This problem has been solved by M. Beck [31]. He found that there still 
exists a genuine buckling load, given by (19.1) with k = 2-031. It is approxi- 
mately eight times as large as the axial buckling load. 


‘ 


24. Pulsating Compression 


If the rod is acted upon by a pulsating load 
P=Q+ Scoswi, (24.1) 


the problem is instationary. According to Theorem 18, the kinetic method 
alone is legitimate; besides, the question of physical imperfections deserves 
attention. 

Case 5 has been treated by E. METTLER [32]. The most general motion of 
the rod consists in an infinity of normal vibrations with amplitudes 4,, satis- 
fying the Mathieu equations 


Ant (bn — En cost) A,=0, (2 = re ey (24.2) 


where derivatives with respect to t = wt are denoted by primes. The para- 
meters of (24.2) are 


A> ) rhe! ) En } 
where the quantities Q,, 2, represent the static buckling loads (of all orders) 
and the natural frequencies respectively of the unloaded rod. 

Let the points P,(6,, €,) in Figure 12.2, for given values of Q and S, be 
represented by polar coordinates r,,, %,. Since Q,~%, their arguments #, 
form a rapidly decreasing sequence. At the same time, since A,, ~ 2, their radii a 
increase so fast that e,~mn*. By the straight line «= 6, the upper half 
plane is divided into two regions of which, on the whole, I is stable while IT 
is unstable. Provided Q + S < Q,, all points P,, according to (24.3), lie in the 
region I. Thus, a rod, as a rule, is stable if the maximum of the pulsating load 
does not exceed the static buckling load. If Q +S > Q,, at least P, lies in the 
region II. Moreover, for Q->Q,, the argument #,->2/2, and for OG 
3, > 2/2. Thus, under a pulsating load the mean value of which is near or in 
excess of the static buckling load, the rod, as a rule, is unstable. 

In order to discuss the most important exceptions from these rules, let us 
restrict ourselves to loads deviating but little from their mean values (S<Q). 
Here, it follows from (24.3) that, apart from the cases Q~Q,, the points P, lie 
near the d-axis. Thus, mean loads Q > Q, are critical. If Q<(Q,, however, 
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instabilities are to be expected for 6, = 1/4, 1, 9/4, ... , i. e. whenever 
yp anerad ee. 
= | Laie sm cyte 


These conditions are met by arbitrarily small values of Q whenever 


= a 
2 Migs Digs Ty 8 (24.4) 


Owing to damping and to the fact that the unstable regions near the d-axis, 
except in the vicinity of 6 = 1/4, are extremely narrow, the most important 
solution of (24.4) is wm =24,. This instability is often referred to as subhar- 
monic resonance. 

In Cases 1 to 4, the problem is more complicated [33]. The influence of 
physical imperfections has been examined by E. METTLER [34]. 


-V. Buckling by Torsion 


A. G. GREENHILL [35] has shown that straight rods and shafts subjected 
to torsion are capable of buckling. In order to discuss their instabilities, we 
shall confine ourselves to shafts having two equal flexural rigidities (8 = «). 
A special case of a prismatic rod with different flexural rigidities has been 
treated by R. GRAMMEL [36]. 


25. Axial Torsion 


Figure 25.1 shows a shaft acted upon by a constant couple of moment M. 
Provided that its vector remains axial during deflection, it shall be called an 
axial moment. Such a moment is (Example 2. 5) circulatory unless 
the end of the shaft on which it acts is held in alignment. Thus, 
Cases 1 and 2 of Figure 19.2 indeed are purely nongyroscopic. 
They have, according to theorems 14, a genuine buckling load M, 
which can be obtained by static methods. However, Cases 3, 4, 
and 5 are circulatory. Their critical moments, according to Theo- 
rem 17, must be calculated by means of the kinetic method. 

According to A. G. GREENHILL [35], the buckling moment 
has the form 


Figure 25.1 M - od (25 1) 
= % ; 
Shaft acted - a be? 
upon by a ; soa F ; E , 
couple. k being anumerical factor depending on the constraints and given 


by Table 25. It has been calculated by A. G. GREENHILL [35], 
the author [5] and A. Tr6scu [6]. In Cases 1 and 2, k defines a genuine buckling 
load. In the other cases, it merely marks the transition of the shaft from a 


174 HANs ZIEGLER ZAMP 


stable to an unstable state. It is probable but has not been rigorously proved 
so far that in the latter cases M, is a genuine buckling load. 


Table 25 
Axtal buckling moments 
| Case | 1 | 2 3 | 4 | 5 | 
2-861 


| Factor k 


ih) 
Oo 
o 
iS) 


The stability of rods both under compression and axial torsion has been 
discussed by A. G. GREENHILL [35], R. GRAMMEL (see [4], p. 545), and A. TROscH 
[6]; for a full account of all buckling cases see [6]. 


26. Conservative Torsion 


Since actually, for arbitrarily small values of M, instabilities never have 
been observed, the results obtained in 25 for Cases 3 and 4 disagree with the 
facts. Thus, the question arises whether it be justifiable to assume that the 
vector of the moment M remain axial during deflection. One might as well 
assume that it remain tangential to the deflection curve; however, an attempt 
by J. Morris [7] based on this assumption has shown that (at least in Case 3) 
the instability does not disappear. In any event, it is evident that any a priori 
assumption concerning the direction of the moment vector is arbitrary and 
that, in order to obtain more reliable results, the forces contributing to the 
moment M must be examined [9]. 

If the moment acting on the section z, of a shaft is produced by two constant 
forces having fixed lines of action (Figure 26.1), the moment vector 


M = M(0, 0, 1) (26.1) 


is axial. This case, however, is extremely rare; perhaps it is approximately 
realized in Pelton turbines. In most cases, however, any deflection of the shaft 
is accompanied by a displacement of the lines of action of the forces and, in 
consequence, by an inclination of the moment vector, though the force vectors 
themselves remain constant. 

In Figure 26.1, for instance, the lines of action of the constant forces may 
be free to follow an inclination of the disc. This type of loading is easily obtained 
by means of threads; besides, it is realized (in connection with an axial thrust) 
in axial turbines having two jets. As soon as the shaft is deflected, the tangen- 


tial vector €(x}, yj, 1) at z, is inclined. So is the moment vector M (Figure 25.1) 
which is given by 


’ 


M = M(x,,0,1), (26.2) 
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provided the forces are parallel to the y-axis. A moment vector of this type has 
been called a quasitangential moment. . 

If two couples of this type are combined as shown in Figure 26.2, the 
moment vector is 

fll | hola: 

Such a semitangential moment is independent of the axes x, y and appears 
whenever a uniform distribution of more than two circumferential forces are 
applied to the disc, particularly if this distribution is continuous as in most 
axial turbines. 


P 
ea y 
Figure 26.1 Figure 26.2 Figure 26.3 Figure 26.4 
Axial or quasi- Semitangential Pseudotangential Two quasitangential 
tangential moment. moment. moment. moments. 


In Figure 26.3 the couple is applied by means of a cross bar fixed at right 
angles on the shaft. This type of loading, too, is easily realized by threads or 
in an electric or magnetic field. If My = 2 Pa denotes the original moment and 
@, the angle of twist just before buckling sets in, the moment responsible for 
buckling is M = M,cos#%,, while obviously 


059% <5. (26.4) 


During deflection, 
M = M(x, + y, tand,, 0,1), (26.5) 


provided that the cross bar in the unloaded state and the loads are parallel to 
the axes x, y respectively. Equation (26.5) defines a pseudotangential moment. 

If the shaft is loaded by means of a Cardan link, M is of the quasitangential 
type. In case 5 (Figure 19.2) the buckling load, provided that the moments 
at both ends are quasitangential, depends on the angle 6 (Figure 26. 4) through 
which one of these couples is rotated with respect to the other one. If the direc- 
tion of the forces P is fixed, so is 6. In the case of two Cardan links, however, r) 
depends on the elastic twist which itself is a function of M; this problem, 
which is more complicated, has been solved by W. T. KOITER [10]. 

Since the moments (26.2), (26.3) and (26.5) are conservative [9], we have 

Theorem 26. Semi-, quasi- and pseudotangential moments are purely non- 


gyroscopic. 
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27. Semitangential Torsion 


According to theorems 26 and 14, a shaft loaded by a semitangential mo- 
ment has a genuine buckling moment which can be obtained by static methods. 

The semitangential buckling moments, calculated [9] for the cases of Fi- 
gure 19.2, are given by (25.1); the factor k being taken from Table 27. Since, 
in Cases 1 and 2, there is no difference between the semitangential and the axial 
moment, the corresponding values in Tables 27 and 25 are identical. 


Table 27 
Semitangential buckling moments 
Case | 1 2 3 4 | e | 
Factor k 2-861 2 


1 | 2-168 | 1-564 


Note that the values of Table 27 satisfy Theorem 20a. 


28. Quasitangential Torsion 


Quasitangential moments are treated in the same way as those of 27. 
The quasitangential buckling moments, calculated [9] for the cases of 
Figure 19.2, are given by (25.1), the factor k being taken from Table 28. They 


Table 28 
Quasitangential buckling moments 


2-861 


ho 


| Factor k 0-5 


differ from those of Tables 25 and 27 in Cases 3 to 5 and again are in accord 
with Theorem 20a. 

In Case 5, provided both external couples be quasitangential, the buckling 
moment depends on the angle 6 (Figure 26.4). Table 28 gives the values for 
6 = Oand 6 = a/2 only. 


29. Pseudotangential Torsion 


Pseudotangential moments are treated in a similar manner [9]. In all of the 
cases of Figure 19.2, however, the buckling moment would produce an angle 
of twist beyond the interval (26.4). Thus, the shafts of Figure 19.2 do not 
buckle under pseudotangential moments acting at their ends. 

The stability of shafts both under compression and conservative torsion has 
been examined by M. Beck [37]; the results will soon be published. 
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30. Pulsating Torsion 
If a shaft is acted upon by a pulsating moment 
M=Q2+Scoswi, (30.1) 


the problem is instationary. A kinetic investigation shows that the shaft 
executes torsional vibrations which, in the case of resonance, may be dangerous. 
Besides, the pulsating moment, like the pulsating load of 24, is capable of 
exciting lateral vibrations. 

In Case 5 and under the assumption of an axial moment, they have been 
investigated by W. KELLER [38]. 


VI. Critical Angular Velocities 


For an observer at rest the critical speeds of a shaft are resonance pheno- 
mena [11]. For an observer rotating with the shaft, they are instabilities of a 
gyroscopic system. In a few simple cases it is possible to obtain the critical 
angular velocities by static methods which, thus far, have been used almost 
without exception. According to Theorem 15a, however, this way of approach 
is unsatisfactory; in numerous problems it does not yield all of the critical 
speeds. 


31. Critical Speeds of an Unloaded Shaft 


Figure 31.1 shows a shaft rotating with the angular velocity w. It is 
equipped with a single rotor of mass m, which is assumed to be statically and 
kinetically balanced. In order to simplify the problem, let us neglect the mass 
of the shaft, its compression and twist, the weight of the rotor and its deforma- 


Case / 2 a 
t t 
AE o\< 
Ye oe) 
Gonslrainls| 7 

Ye Up 

Ae dE > 

di ig fs 
eaiowtere ae Wa Wa 

Figure 31.1 Figure 31.2 


Shaft equipped with a rotor. Different cases of critical speeds. 


tion as well as frictional forces and gyroscopic moments. Then the flexural 
vibrations, referred to a coordinate frame rotating with @, are determined by 


m (%—20)¥,—0?%)=—-F,, mi, +20%— 0? y;)=—h, (31.1) 


ZAMP 1V/12 
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where F,, F, may be interpreted as the forces producing the static displace- 
ments Xj, Vy, of section 2. 

Provided that the flexural rigidities «, B of the shaft are equal, the deflec- 
tion a at z, due to a unit force acting at z, is parallel to this force. It has the 
character of an influence number and depends on the constraints of the shaft. 
Introducing this influence number which, in three simple cases is given by 
Figure 31.2, into (31.1), we obtain 
a ee (ae ~o)%4=0, t2o%+ (4 -#%)%=0. (31.2) 
Hence, the problem is purely gyroscopic. Its differential equations are those of 
Example 8.1, provided the particular case c, = c, = 1/a is considered. Thus, 
the shaft is stable for any value of w. 

In order to explain the existence of a critical speed, physical imperfections 
such as static or kinetic unbalances must be taken into account. They can be 
represented by a few small additional masses located at arbitrary points of the 
rotor, giving rise to constant inertia forces. These forces, according to Theorem 
lla, produce a latent instability whenever a nontrivial equilibrium position 
exists. Thus (and in accordance with Example 11.2), the angular velocity 


(31.3) 


1 
0) = —— 
he Ym a 
is critical (if negative values of w are neglected). 
32. The Influence of Gyroscopic Moments 


When the gyroscopic moments of the rotor are taken into account, the 
problem still is purely gyroscopic. 


= —05 0 O05 / 
: Figure 32.1 
Influence of a gyroscopic moment on critical speed. 
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It has been treated by A. Stopora [39], R. GRAMMEL (see [40], vol. II, 
p. 13) and the author [11]. Here again, the critical angular velocities represent 
latent instabilities. Provided the rotor is flattened (its axial moment of inertia 
C exceeding the equatorial moment A), there appears a single critical angular 
velocity w, > w,. If the rotor is elongated (C < A), there exists a pair of them 
(one greater and the other smaller than a). 

In Figure 32.1, @,/@ is plotted versus k?/l?, the quantity k? = (C — A)/m 
being positive or negative according to whether the rotor is flattened or elon- 
_ gated. The curves are obtained by using the influence numbers of a cantilever 
shaft; thus, they are valid in Case 1 of Figure 31. 2. 


33. Shafts Having Two Unequal Flexural Rigidities 


Provided the flexural rigidities «, B of the shaft are unequal, the elastic 
force (F,, F,) introduced in 31 generally does not have the direction of the 
deflection (%,, y,). If the influence number a,, (7, k = x, y) represents the 
deflection at z,, measured in the direction 7 and due to an unit force parallel 
to the axis k, we have 


GeO Ge = Opn O Ayye 0 (33.1) 
and 
1 1 
i, = aie 4, i= on, yi: (33. 2) 


Neglecting the gyroscopic moments again, we obtain, by introduction of 
(33.2) into (31.1), 


#,- 205, + (2— - 0°) 4 = 0), ji + 20 % + (2-0!) y= 0. 


M Ay» m ayy 33 3) 


Thus, the problem still is purely gyroscopic. Its differential equations are those 
of Example 8.1; hence, the angular velocities 


ae teres ee Se Se ae a) 1634) 
VM ayy 


= Ym ay, 
are critical. 

This result has been established by R. GRAMMEL (see [4], p. 786). According 
to (33.3), both limits of the critical interval (33.4) represent latent instabilities 
and might have been obtained by static methods. According to Theorem 15a, 
however, it cannot be proved by static means that the whole interval 1s 
critical. Hence, the principle of equivalence which, likewise, is due to R. GRAM- 
MEL (see [4], p. 782), stating that the critical values of m are those for which, in 
the absence of unbalances, nontrivial equilibrium positions exist, is not gener- 


ally valid. It must be replaced by 
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Theorem 33. In a purely gyroscopic system, the angular velocities w, < Wy, < ... 
admitting nontrivial equilibrium positions (in the rotating coordinate frame), are 
critical. Any other angular velocity w > w, may or may not be critical. 


34. The Influence of Compression 


The shafts of Figure 31.2, loaded by a constant axial thrust P, give rise 
to the five cases illustrated in Figure 34.1. Since P is conservative, the problem 


| 
: 
| 


Figure 34.1 
Loaded shafts. 


still is purely gyroscopic. Provided that « = f, the differential equations of © 
motion are (31.2); the critical angular velocity therefore is | 


(34.1) 


1 
oO, =——, 
7s Ym d 
the influence number @ denoting the deflection at z, due toa horizontal unit force 
at z, and the vertical thrust. 

The evaluation of a shows that, for loads small compared with the buckling 
load, the ratio of m, and qs» is given by the formula of R. MELAN [41], 


Oy, _ 0 P|? 
ot tae (34.2) 


M9 
the factor k being taken from Table 34 (see also [4], p. 800). 


Table 34 
Influence of a thrust on critical speed 


35. The Influence of Torsion 


If the shafts of Figure 34.1 are acted upon by a torque in place of the 
thrust, the differential equations (31.1) remain valid. Due to the torque, however 
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(and in spite of the assumption « = f), the force (F,, F,) does not necessarily 
have the direction of the deflection (x,, y,). Making use of the influence numbers 
defined in 33, we have 


fom Gate day Fes = aie ay, F (35.1) 
The determinant 
eer A Any Cag lise (3562) 
is positive (see [4], p. 89) ; hence 


1 u 
Bit a A (yy %— Agy 1) ’ EF a A eae xy i Ay a 44) 5 (35.3) 


Introducing (35.3) into (31.1), we obtain the differential equations 


m A 


- E Ga be 
%—2wyy+ (= = wo) mg at V0, | 


. (35.4) 
Py + 2 ty — ea + (08 — wt) = 0, | 
and the characteristic equation 
a (2 w? + “ee Svs) A2 4. 2 ew Sve gy 4 | 
\ M Z mA 
é (35.5) 
4 Gen T Uy 2 " 
ae mA “a m? A : 


If all of the couples acting on the shaft are axial, the rotating system 
(except in Case 2) is not only gyroscopic but also circulatory. This implies, ac- 
cording to 10, that a,, + dgy. 

In Case 1 the shaft at rest, according to 25, is unstable under arbitrarily 
small values of M. It has been proved in [5] that it remains unstable for any 
value of w; thus, any angular velocity is critical. A. TROscH [6] has shown that 
the same result applies in Cases 4 and 5, while in Cases 2 and 3 the critical 
angular velocity (for moments small compared with the buckling moment) is 
given by the formula of R. GRAMMEL (see [4], p. 800), 

Moya (2), (35.6) 


Wo a 
the factor k being taken from Table 35a. 


Table 35a 
Influence of an axial moment on critical speed 


Case 


iS) 
Oo 


Factor k 40 oR 
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In most actual cases, however, the moments acting on the shaft are either 
semi- or quasitangential. In Case 2, this is irrelevant. In Cases 1 and 3, the ~ 
moment acting on the disc must be specified; in Cases 4 and 5 those acting on 
the ends must also be specified. 

If the moment acting on the disc is semitangential (as is the case in turbines 
having more than two jets), the problem, in Cases 1 to 3, is purely gyroscopic ; 
hence 


Gnu 05 A Cee (35.7) 

Moreover, since the semitangential moment (26.3) is independent of the coor- 
dinate frame, dye = Ayy and aye = —az,y; thus 

Ong = hyy= 470, Bey =4,=9- (35.8) 


It follows that (35.5) is quadratic in A?. The discriminant of (35.5) is 


Dae ert tn gt y (Aee + yy)? = 44 (35.9) 


m* A* 


and the roots 42 are negative or zero. Hence, the only instabilities are of the 
latent type. They are given by (34.1), @ denoting the deflection at z, due toa 
horizontal unit force at z, and to the moment M. 

The evaluation of @ shows that, in Cases 1 to 3, formula (35. 6), together 
with the factors k taken from Table 35b, is valid [11]. 


Table 35b 
Influence of a semitangential moment on critical speed 


Case i 


Factor k ates : eS 


In Cases 4 and 5, the moments acting at the ends of the shaft still may be 
of various types. If applied by means of Cardan links, they are quasitangential, 
i. e. of the form (26.2). The problem, which is being treated by CH. WEHRLI[42], 
is purely gyroscopic. Since the relations (35.7) still hold, the discriminant 
(35.9) is positive; hence, the roots A? of (35.5) are real. They are both negative 
if and only if the constant term in (35.5) is positive. Thus, an entire interval 


O, Sw SW, (35.10) 


is critical, the limits of which follow from 


1 
2mA 


(dex + Gyy FV (dee — 4yy)? + 4 a2, )- (35.11) 


9 
M7, = 
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In Case 5, Ca. WEHRLI [42] obtained the critical interval 


aM : 
Tee ee Sept! 1a 
64 leg = Wo 64 “ I 


(35.12) 


Case 4 is more complicated, since the critical interval depends on the angle 6 
(Figure 26.4) which, by itself, is a function of M. 
_ Ifthe moment acting on the disc is guasitangential (as is the case in turbines 
having two jets), (26.2) holds in a coordinate system at rest. Therefore, the 
problem is best treated in a fixed coordinate frame. 

In Cases 1 to 3, this treatment [11] shows that there are two critical angular 
velocities. In Case 1, they are given by 

uO ta eee. (35.13) 


Wo 16 a 


In Cases 2 and 3, they coincide for small values of M, being given by (35.6) in 
connection with Table 35c. 


Table 35c 


Influence of a quasitangential moment on critical speed 


Case | iz | 3 
Becton B 3 13 
ee 640 7680 


Case 4 does not present any new aspects, since the disc is not loaded. In 
_Case 5, the assumption of a cardan link at one end implies that at the disc (26. 2) 
holds with respect to a fixed coordinate frame, at the end, however, with 
respect to a rotating one. As a consequence, there exists no coordinate frame 
in which the influence numbers are independent of the time. The problem 
which, in this case, is instationary, is still unsolved. 


36. Conclusion 


The problems treated in 19 to 35 are taken from a comparatively na. 
field. Yet they clearly show that the range of validity of static methods is 
more limited than is usually supposed. Besides, it seems probable that. this 
result will have some bearing on other problems treated so far by static 


methods. 
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Zusammenfassung 


Elastische Systeme k6nnen im wesentlichen in konservative (und zwar 1. nicht- 
gyroskopische und 2. gyroskopische) sowie nichtkonservative (3. dissipative, 
4, zirkulatorische und 5. instationare) Systeme eingeteilt werden. Eine Analyse 
der fiinf Kategorien zeigt, dass die bei Stabilitatsuntersuchungen tiblichen sta- 
tischen Verfahren nur in der 1. und in beschranktem Masse in der 3. Kategorie 
zulassig sind, wahrend die iibrigen Systeme kinetisch behandelt werden miissen. 


(Received: January 7, 1953.) 


Beitrag zur Kenntnis der Eigenschwingung einer idealen 


Fliissigkeit in kommunizierenden Rohren 


Von Hsten-Cu1H Liv, Schanghai}) 


1. Einleitung und Geschichtliches 


Isaac Newton?) behandelte die Schwingungen von Wasser in einer U-for- 
mig gebogenen Réhre mit senkrechten Schenkeln und konstantem Querschnitt. 
Wenn L die Lange der Wassersaule ist, so betragt die Zeit einer vollen Schwin- 


gung Hs lee (1) 


JoHANN BERNOULLI*) zeigte, dass, wenn die Schenkel des Rohres die Nei- 
gungen a, und % gegeniiber der Waagrechten haben, die Schwingungsdauer 


re 2, 
T= an|/- (sin a, + Sin a) © 
ist. 


DANIEL BERNOULLI) 5) untersuchte die Schwingungen von Wasser in ver- 
bundenen Gefiissen mit senkrechten Schenkeln. Es werde angenommen, dass 


1) Technische Hochschule Tsinan. 2 
2) Isaac Newron, Philosophiae Naturalis Principia Mathematica (London 1687). 
3) JoHANN BERNOULLI, Commentarti Academiae Scientiarum Petropolitanae (1727). 
4) H. pe LAGRENE, Cours de Navigation intériewre (Paris 1869). 

5) DANIELIS BERNOULLI, Hydrodynamica (Argentorati 1738). 
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die Endriume Zylinder bilden oder, bei beliebiger Gefassform, dass die Schwin- 
gungen in derart engen Grenzen bleiben, dass man /, und Ff als konstant ~ 
betrachten darf (wobei man unter /, und £, die Mittelwerte der Querschnitte 
im Bereich der Schwingung versteht). Ferner sei 


= — = const 


aes L 
Be Aer 


gesetzt (wobei J, der Querschnitt des Rohres an einer beliebigen Stelle und ds 
ein Linienelement langs der mittleren Stromlinien bedeuten), namlich bei zylin- 
drischen Endréiumen /, = /, oder nahezu konstant, namlich bei ungleichen 
zylindrischen Endraéumen. Auch sei jeder Zylinderquerschnitt viel grésser als 
jener des Verbindungsrohres oder auch bei beliebiger Gefassform die Schwin- 
gung recht klein. Unter diesen Annahmen erfordert die Vollschwingung die Zeit 


hee 
Brae poke Ba 
P= 22\ rey a 


Das Ziel der vorliegenden Arbeit ist die Untersuchung der Eigenschwingung 
einer idealen Fliissigkeit in kommunizierenden Réhren von veranderlichem 
QOuerschnitt. 


2. Ableitung der Bewegungs¢gleichung 


Nach der erweiterten Bernoullischen Gleichung fiir nichtstationare Stré- 
mung ist fiir die beiden Spiegelquerschnitte wegen ~, = f, = Atmospharen- 
druck 


2 
ve 


Se 
; Us ’ Lie 
— $s, sing, = So SiIN& + — | ds. 
2¢ 1 hea + Sy SIN &, + mall ~ ds (4) 


Fig. 1 
Festlegung der Koordinaten fiir das Schwingungssystem. 
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In dieser Gleichung und in Figur 1 bezeichnen v, und v, die Spiegelgeschwindig- 
keiten, v die Geschwindigkeit an einer beliebigen Stelle s zwischen den beiden 
Spiegeln, s, sing, und s, sing, die entsprechenden lotrechten Komponenten der 
Verschiebungen der Spiegelflachen gegen die Gleichgewichtslage, und ds ein 
in Strémungsrichtung gemessenes Langenelement. Wegen v, i =v, und 
5,24 = Ss, Fy folgt aus (4) 


we F, \2 : Sie’: Te! ov 
ea [1 — a = (sina, + sin as] = . ii Te ds 
und unter Beachtung von v, = ds,/dt und vF = v, wird 


[1 — ()"] esi —2¢ (sin 4, + oe sin as) eS ss / ae ds 


yy) \ F, eI I da 
oder mit . 
F, ‘Ce ies a ean tes ops 
= 5" SS = = S6S=L und 2A=72— 1 
jee hae fp | z 
Sh Ss 
2 ib ds, \2 : : 
fe — ae ( a) + g (sina, + Sing) s,=0. (5) 


a, Umformung der Bewegungs¢gleichung fiir den Fall F, = const, 
F, = const 


Die beiden Schenkel der Rohre (Figur 2) mégen verschiedene Neigungen 
besitzen, also a, + %. Jetzt gilt 


lov tap peeane é 
=a — 1 =const. (6) 


Spezialfall des Rohres; 0 — 0’ die Gleichgewichtslage; F =; = const bei Stelle 1; 
F = F, = ‘const bel Stelle 2. 


Fiir die reduzierte Lange L wird zweckmissig geschrieben 


s 


. 


nds + | nds; 


pi = [yas = | nas — 
Sy 0 J 


S20 Sy Ss 
~ 


ae = 7 ds = const, [a dS = 81, fu ds = 92 Sy, 


0 0 S20 
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hierin sind 


also L = L,+As,, und damit lautet die Bewegungsgleichung: 


\2 


(Ly +As,) J +g (sinay + mp sina) 5 = 0. (7) 


Ae Nag [ids 
dt? + 2 A | dt 


4. Integration der Bewegungs¢gleichung 


Der Wert 7, ist je nach der Rohrgestaltung sehr verschieden, und demnach 
kann auch die Lésung dieser Differentialgleichung sehr mannigfaltig sein. Die 
Untersuchung soll hier nur auf den Spezialfall A = const, & = const beschrankt 
werden, da man sonst auf grosse mathematische Schwierigkeiten stésst. Es 
wird also vorausgesetzt, dass die beiden Enden der Flissigkeitssaule sich stets 
in Rohrstiicken von konstantem Querschnitt bewegen. ; 

Es seien 2 $;max (Figur 2) die Lange des Stiickes 1 mit A = const und 
2 Sy max die Lange des Stiickes 2 mit A= const. Das Rohr kann durch die 
Ziffern I, II, III und IV in drei Abschnitte unterteilt werden, und zwar 


IV HI 
” a 
L,= / 7B 4 = Syman + / Sr 4S + he Seman 


- 


I Il 


Nach der oben erwahnten Einschrankung sollen die Enden der Flissigkeits- 
sdule wahrend der Bewegung sich niemals in dem Abschnitt II-III befinden. _ 
Infolgedessen muss sein: 


Sig <f “Sree 4 (8) 


Nach Multiplikation von (7) mit s, lasst sie sich in die folgende Form bringen 
d he ‘o 1 *) * s 9 
7 ( = si + z As s+ 4 (Sin a + M2 SiN a) si) a= () (9) 


oder nach Integration 


ds, \2 5 : 
(Ly +A 5;,) (3) + g (sina, + 7 sina) s? = A , (10) 
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wobei 


[4 = g (SING, + Hp SING) Sip (11) 


durch die Anfangsbedingung ¢ = 0, s, = 0 und s, = 5,9 bestimmt ist. s,) bedeutet 
die maximale Anfangsamplitude. Der Ausdruck (10) lasst sich schreiben 


- ——— 2 2 
eee = Vg(sing, + No SiN oy) = : (12) 
Gleichung (12) besagt, dass die Maximalausschlage auf beiden Seiten der 
Gleichgewichtslage gleich sind. Auffallig ist das Auftreten von s, im Nenner 
unter der Wurzel. Daraus lasst sich schliessen, dass die absoluten Werte der 
Geschwindigkeit, bezogen auf gleiche Abstande zu beiden Seiten der Gleich- 
gewichtslage nicht gleich sind. Wenn man s, gegen s, auftragt, so bekommt 
man ellipsenartige Kurven (Figur 3), deren Maximalwert jedoch nicht auf der 
Achse s, = 0 liegt. Nach welcher Seite und um welchen Betrag dieser Maximal- 
wert von der Mittellage abweicht, wird durch den Wert 7 bestimmt. 
Bei dem Ausschlag 


(13) 


Fig. 3 
Geschwindigkeiten und Beschleunigungen in Abhangigkeit vom Wege fiir 1) = 7. 
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hat die Geschwindigkeit ihren Maximalwert, wobei das positive Zeichen von 
der Wurzel fiir 7, > 1 und das negative Zeichen fiir y, < 1 zu setzen ist. 
Man erhalt noch 
ds, \2 z 2 : Per Ss 
(+) = Ga & (Sina + Np SiN %) (5 +VL3—- so) , (14) 

wobei das negative Zeichen fiir 7,>1, das positive Zeichen fiir 7, < 1 gilt. 
Ferner ist die Beschleunigung 

as; 


: . 7 
or ear (Sina, + 72 SiN%,) — 


Sip + 2155, + 457 
2 (Lo + A5,)? 


(15) 


beziiglich der Gleichgewichtslage unsymmetrisch (Figuren 4 und 5). Bei der 
Gleichgewichtslage hat sie den Wert 


d-s, 
( dt? 


Jno ~€ Con + sine) > (7A) (16) 


Fig. 4 


Wege, Geschwindigkeiten und Beschleunigungen in Abhangigkeit von der Zeit fiir No = 7. 


| 


: 


4) URS ; ; F ; ae ee oe 
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Fig. 5 


Die maximale Beschleunigung in Abhangigkeit von den gréssten Auslenkungen. 


Die maximale Beschleunigung hat den Ausdruck 


( os = § (SiN % + Mg SIN %) Smq~ (17) 
at | max ae ati (73 i) Smax 


Nach Trennung der Veranderlichen und Ausfiihrung der Integration nach 
der Zeit in (12) ergibt sich 
/L + FAS, Bi 


Sy 
| ie (18) 
Ve ( sin a, ae Na ‘sin Op) Sto — S} 


ft, = 


Setzt man s, = sj), cosy in Gleichung (18) ein, so erhalt man 


Po a 
Rei ) Jap 
tale = 2 (sin a, + 7 Sin a) "| Asie + cosy dg . (19) 
Pi . 


Nach weiteren Umformungen lasst sich die Gleichung (19) in der Form schreiben 


YW. 


— Va Lo +A 51 i V1— k sin? y dy; (20) 


sin % + 1 SiN &) 
Ya 
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darin bedeuten : 


Die rechte Seite von (20) ist ein elliptisches Integral zweiter Gattung. Der 
Ausschlag sei willkiirlich als positiv festgelegt, wenn die reduzierte Lange mit 
wachsendem Ausschlag zunimmt, und umgekehrt. Die Schwingungszeit von 
Null bis s,) und wieder nach Null zuriick wird 


[Lorik sy ds, 


= 2 
RESET nL fV Sio— Si 


le - 


Mit Ss; = Sy COS@P: 


eee - IE 
Ly. 10 0 
ards V3 (SIN % -+ Jz SIN % ) / | A S10 a chet 
0 


oder mit p = 2 yp: 


= ot A Sio = aT RC Ke 
Toe = 4 rat [vie sin? y dy , 


oder schliesslich mit der iiblichen Darstellung fiir elliptische Integrale zweiter 


Gattung 
Lg tA So Dkeee x 
fae 3 (Sin %, + M9 SiN oy) B(/7 yy Wee ey 2 (21) 


Die Zeit fiir die Halbschwingung in der entgegengesetzten Richtung von Null 
bis —sy9 und wieder nach Null zuriick ist 


=. NV ot4s 2 
gitin Cee si—st 


le = 


mit Ss; = Sy) Cos@: 
Pee AS 
T,,=2 i Jee 
a Vz g (Sin a + M9 SiN ay) | Vis + cosp dp 


oder mit p = 2 p 


} he n/2 
Wi oenstar feel 
Tw» = (Sin % + my Sin a) | V1 — k* sin*y dp 
n/4 
oder schliesslich 


Loraty eh RIC 
ap = //, 10 10 4 
a — 4 g (Sin % + Mp Sin ag) |E E ( "Ly Ase 4 “ale (22) 
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Die Zeit einer Vollschwingung ist 


Lot S1o 
SIN & + 15 SIN &) 


T* = T+ Tw = 4/5 E, (23) 
6 

worin £ ein vollstandiges, elliptisches Integral ziweiter Gattung ist. T* ist ab- 
sichtlich eingefithrt, um den Fall hervorzuheben, bei dem die reduzierte Lange 
wahrend der Schwingung nicht konstant bleibt. Durch Reihenentwicklung und 
nach gliedweiser Integration erhalt man mit 


Te 2 |/ les +4 S19 1 = 1 WE ZASio) \t > 
<a g (Sina, + 72 Sin a) 27 ea Vee Gen rel - (24) 


Es lassen sich noch die folgenden Verhdltnisse ableiten: 


T* E Te BE (h, nA) T* E 
Pet ee Fee Ty Stee 
Da nach der Eigenschaft der elliptischen Integrale stets 
E(k, 2) <E<2E(k, =| 
gilt, so kann man noch die nachstehende Tatsache folgern 
Tow < Tin (26) 


(giiltig fiir 7, +1). Fiir Schwingungen der Fliissigkeit in kommunizierenden 
Rohren, bei denen 7» +1 ist, sind die zwei Schwingungszeiten T,,, und Ty, 
in beiden Schwingungseinrichtungen verschieden (Figur 6), und zwar ist die 
_ Schwingungszeit in der Richtung, in der die reduzierte Lange mit Zunahme der 
Auslenkung wichst, grésser als die in der entgegengesetzten Richtung. 

Unter Beibehaltung von L, setzt man in (24) 7, = 1, dann ergibt sich 


Lo 
Beet ls (sina, + sin a@,:) ~ 


Das stimmt iiberein mit Gleichung (2), die von JOHANN BERNOULLI bereits 
abgeleitet worden ist. Wenn man weiter a = % = 2/2 setzt, dann wird wie in 


Gleichung (1) tual 
Lo 
I ee Opens ee . 


Bei gleicher reduzierter Lange Ly erreicht 7* ein Minimum, wenn die beiden 
Schenkel des kommunizierenden Rohres lotrecht stehen. Fiir ein senkrechtes 


ZAMP 1V/13 
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Fig. 6 


T*, T,,, und T., in Abhangigkeit von 7, mit der Anfangslenkung s,, als Parameter. 
ab 5 2 § § S10 


, zu 


U-Rohr mit ungleichen Schenkelquerschnitten /, = const und & = const ist 
die Schwingungsdauer als Funktion von /j, /5 und 24 (Figur 2) 


si | Lot (n5 ; 1) wo f _> um 2(n% = 1) m9  @ 
) n=1 ‘ ; 


Tyy = 4 =e) 20 E ~ Ea eS (29) 
Fir sehr kleine Ausschlage kann man die folgende gute Annaherung erhalten 


|/ Fy Ly 


THs eA) } 
T* ww 20 |) oe (30) 


(giiltig fiir sehr kleines x,)). Setzt man 


= Pa ee 
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so erhalt man 


aA | ae ea 


was wiederum mit Gleichung (3), dem Resultat DANIEL BERNOULLIS, tiberein- 
stimmt. 


5. Vergleich zwischen T und 7* 


Wir definieren T als die Schwingungsdauer einer F lussigkeitssaule in einem 
kommunizierenden Rohr von konstantem Querschnitt. Die Lange dieser Fliis- 
sigkeitssdule sei Ly. Dieses System soll mit einem zweiten kommunizierenden 
Rohre verglichen werden, das jedoch ungleiche Schenkelquerschnitte aufweist. 
Die Lange der Wassersdule sei auch in diesem zweiten Fall — in der Gleich- 
gewichtslage — gleich Ly. Die Neigungen der Schenkel seien fiir beide Falle ent- 
sprechend gleich. 7 soll als ein Vergleichsmass fiir 7* dienen. Von den friiheren 
Resultaten lasst sich bestatigen, dass 


ie 


Si IG SIN & + Np SIN o& Tt 
AU -|/, Sin % + SIN Oy _ 
Sin a + 7» SIN om 


Da 7 > 1 ist, gewinnt man offenbar die Beziehung 


(a arte Lot (n3—1) sy Sina, + sina, ace 
Sip 0 


Lak (31) 


Bee) 


Die Bedingung 0T*/0s,, = 0 liefert s,) = 0 als Lésung von 
A oe a eet Enel Re i oe 
We ris! Dyes | | DO Ae aoe ; 


Mit sj) = 0 kann man zeigen, dass 


(a Te (43 — 1)® 


Ost es 3 Vg (sin a+ 12 Sin a9) 4 15" 


Zufolge dieser Untersuchung besitzt 7* bei s,) = 0 ein Maximum, woraus man 
schliessen kann, dass stets gilt 


TEST (32) 


(giiltig fiir 7, +1). Die (7*, sy9)-Kurve besitzt bei s,)= 0 eine horizontale 
Tangente und fallt mit wachsendem sj) monoton gegen die syy-Achse ab. 
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6. Eigenschaften von T,, und 7,, 


Aus der Beziehung 
OT, 


“ZU 


(2) a! 
OS19 /s1-0 Vg (sin @, + No Sina) Ly 


ist ersichtlich, dass die (T;,,, S49)-Kurve bei sj) = 0 nur horizontal verlaufen 
kann, wenn 7 = 1, das heisst A, = & ist. Wenn 7, > 1 ist, steigt die Kurve 
mit wachsendem sj, und wenn 7 < 1, fallt sie monoton ab. 

Die Eigenschaften von 7,, kénnen aus der Differenz 7* — T,,, abgeleitet 
werden. 


7. Vergleich zwischen 7/2 und T,, 


Aus (21) erhielt man 


. 
Lo 
(Taulec>0 ahd V; (Sin a + 12 SIN &) : 
und 7 
q 
tT oA f = ne eae 
ae ies g (sina, + sing.) * 


Man kann daraus schliessen, dass 


(Teu)s,-+6 a (33) 


n]y 


Bei zunehmendem sj, ist aber auch 7,,, > T/2 méglich. 


Summary 


The oscillation of an ideal liquid in communication tubes of arbitrary shape 
leads to a non-linear differential equation of second order. The coefficient of the 
second derivative is a linear function of the displacement, and furthermore the 
differential equation is not linear, because a quadratic term in the first derivative 
appears. 

The equation becomes linear if and only if the cross sections at the ends of the 
liquid column are constant but not equal, the periods of oscillation in the two 
directions are different and the difference increases monotone with increasing 
amplitude of the displacement. 


(Eingegangen: 27. November 1952.) 
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Torsion of a Circular Shaft with Diameter Varying 
Periodically along its Length 


By Hajimvu Oxkupo, Sendai, Japan?) 


In this paper, the torsion problem is treated for a shaft of circular section, 
the diameter of which varies periodically along the length of the shaft. An 
approximate solution for the problem is derived on the assumption that the 
dimensions of a groove are small as compared with those of the shaft. In the 
assumption used here, however, the dimensions of a groove are meant to be 
comparatively small, and not meant to be infinitely small as is assumed in the 
work of NEUBER?). . 

We may take the axis of the shaft as the z-axis, and use the polar coordi- 
nates (7, 8) for defining the position of a point in the plane of a cross section. 
Then, the stresses, which do not vanish, are expressed by the equation?) 


G oy G dy 
Ci Ge tph Op aa OUE. a saat el a (1) 
where G is the modulus of rigidity, and the function y satisfies the differential 
equation 
0? y a oh Ory 
ee a Te ee (2) 


As the outer surface of the shaft is free from traction, the function must satisfy 
the condition 
wy = const (3) 
at the boundary. 
As the stresses for the present problem vary periodically with z, we take 
for the function y the expression 


poten DA, I, (=2") cos “=, (4) 


n=1 


where 6 is a constant, J, is a Bessel function, and 2a is the period of the grooves 
along the length of the shaft. Let v denote the displacement normal to the axial 


1) Tohoku University. ’ ; 
2) H. NeuBer, Kerbspannungslehre (Springer, Berlin, 1937), p. 140. — 
3) S. TimosuenKo, Theory of Elasticity (McGraw-Hill, New York, 1934), p. 276. 
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plane. Then we have 
Omen ue Ow 
3 _ = —— 4 a 5 
Oz ( ) Ole (5) 


Substituting the foregoing expression for y into (5), we obtain 


0 v gta ee gaily NZ 
+5 (=) =0+ ay 2 ALG Pee sare. 


Integrating both sides of the equation with respect to z over a period of grooves, 
it follows that 


Sle. rr [=]... = 246. (6) 


The left hand side of equation (6) denotes the difference in the rotations between 
two cross sections z = + a. Consequently # is the mean twist in a period of 
grooves. 

Let us now assume that the dimensions of the axial section of a groove are 
small in comparison with those of the shaft. At the points where + is large — 
compared with a, the function y is given by the expression 


0 ail2y/2 SY 4, 7 N12 
= 4 n ‘in j nara xo 
aera lin lakes ok ’ (2 nz ° COR STs (7) 
n=1 
by virtue of the approximate relation 
ny Het IT We LS 15a 
i ( oo ) oe where 7», = 1—— =— 
a a a(2nr)il2 ? In 8Snayr 


Let us consider another function y*, given by the equation 


yt=prtt eo? SCA Ripal, 8) 
n= 
where R denotes the real part of a complex variable, C,, is areal constant, and 
Pi = log (1+ e™*-*),  o = = faa p Ps eS = f ouat, vee (C= 74242) 
If ry = 1, y, may be expanded in the series 


ee (<2tjeet ss d 
een hae | \ : ee nx(C—1)/a 
Py alg y+ an e : 


except at the points y = 1, z = (2 + 1) a, in which 7 is an integer. We take 
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for m, the expression 


Gi 2 0 7 (tt % 
a / PPA Thee = —nm(C—1)/a 
Pe 2! a? (5 1) AG pd n2 é 5 
m=1 
Ae z mee 2 i ahaa x 
P3 = zr (¢ ext 1s aE & (C 1) ) 7 ( 2 ; ge Nn(S—I)/a ; 
; : n=1 
qe pad RN eck 7a Sy (—1)"41 : 
—— / 455i IGG 2 Brat a —n2(f—-1)/a 
Pa 4! at (¢ —1)*4 2! 6a? (S = 1)? 360 oe, n4 G 2 


Pe Sn Oe SMS ae Combis) 916 el edie. saelie! @ (9) je. k6\ 9's 6 s%1ei'6 se. 0°'erhs, 9: ois) (e et mumewts-'e Te “e wie © eile (oie ee ef ey erie WS) we ce 


To simplify the discussion, let us consider a portion of the shaft bounded 
by two sections z = + a, and denote the value of z in this portion by 2. Then, 
the foregoing expansions for q,, hold for all values 
of z= 2% +2na, consequently R[y,] is a many- 
valued function of z. In the calculations for R[,], 
however, we must take z) for z, since the shearing 
stress T, on the section %) = 0 must vanish from the 
symmetrical condition of the shaft with respect to 2p. 

If the equation 


_ =l+h(a—z2), (@2220) (9) 


represents the outer surface of the shaft, Figure 1, 
and the function »* satisfies the boundary condi- 
tion (3), except at the point y= 1, z =a, it follows that 


S{1t+h(a—2h+{ith (a a} COR O,\= 6, (a> 2 Oe AO) 


where 
R[¢:] = <= (a—2z) + Y (SGN Eisen eee ace 
Rigsl= sir — {h? (a — 2)? — 22} 42 oy (-1)"*2 pra st cos BEE, 
Rls] = sige <= {h8 (a — 2) —~3h(a—2z)z )2}4 2 neh baa 
e > ( Iyer Awe: 1 ala) lh 


C0) 2 0 WKS ei pie 6 4 bie 
aianelisy siete (@ lw. 9: spaharsie ie s @ eer © 
Dialicteiiete wise shoelace 600. sf0 4 ©, 6) ¢ 6 6 lee se, e%e = 
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As ais a small quantity in comparison with unity, omitting the smaller quan- 
tities of the order a*, we have 


{14h(a—2)}? =(1+ ah)? —Zh(l+ah)z + => W(t ah)? 22, 


ty] vo re] un 


{1+h(a—z)} P= (1+ ah)-?? + —h(l+ah)7? 2+ Se h2 (1+ ah)-7? 22, 


Using the foregoing approximations, we rewrite the equation (10) in the form 


Le °] 


Oo + Oz + Ogz2 tee Dy P,, e-M*#(1-218) cog ee wi 3 (11) 
where 
2 aan. feng Ps, SERA Css Py ee Cust 
ig Bm a 
Q,= ane 6h (1+ ah)>? — a h(ltah)-? C+ ae Py or C, 
i ts pa ear Cue tt 
Qn = 42 OWE (1-4 ah) — 3 21 + ah) G+ 2 > - or C5 
Site BEM cae 
ana SEM ae 
afi pues \ Cy en = RUG — A"). 


= m x a 
“ , mM IMZ 
2 = y Ap, COS > t= y A», COS : A 
m=O m=0 
h NMZ ses m 2 
gonah(l—z/a) ... CLS nae TZ 
; COs ae i by m COS a , (12) 
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where 
a a2 
4 =~ Ds , A a ae d ? 
2a i 4 a2 (—1)™ 
ome mea 1)" ie Ce a ee , (m = 1) 
ea ear a hh paces. 
ae I (1 + h?) ; 
h [(—1)"—™ = e-nan (—1)rt+m _ e-nah 
ag ee Sa 
nm 1 n2 f2 +- (n a= m)? n2 2 + (n a m)? | : (m = 1) 


Substituting the series (12) into (11), it follows that 


e+ 3 (ami + a Qe POU Sie cos" =x ()y, (13) 


from which we have a system of linear equations for the determination of the 
unknown coefficients C,, in the form 


Qo 4 Oy 4 Oe +o + D? Ong Pa = 9, 


dy Opt @, Op Fort Dy bay Pa = 0., (14) 
n=1 Pp 


ot ie! “eieuny © lel’a: ‘aise! 6ilie! ecm 6) Ome: 6) (016 LO 46) 10) G, 6) &. 6) OL ei 0) @) «Tener 8: 6 


Taking the finite terms in each series in (14), and solving the simultaneous 
equations, we can find approximate values of C,,. Thus, we can readily find 
the corresponding stresses for x = 1, from (8), as 


& 


. A a) 
FG Tk G YG, R|> ae reel | 
ex (15) 


peer oy C, RII “| | 


Let t, be the shearing stress along the outer curve of the axial section of the 
shaft. Remembering that the shearing stress normal to the bounding curve 
vanishes, we have 


(1 + h?) y1re G (1+ Pn 
= OE Deh area 2 oR Ge]. 08 


where the subscript s denotes the value at the boundary. © 
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If y <1, the functions gy, @»,--. are expanded in the series of different 
forms, except at the points” = 1,z= (2 +1) a, as 


foe) 
(= 1)ers —n2(1—C)/a 
Coie ys n2 é ? 


ca (—1)"* —nx(1—)/a 
P3 rT é , 
=1 


Hence 


p* = = yt ds (—1)"+1 4’, e714 cos se J (17) 


where 


ioe) 

A‘ a enna Cm 

Ae te mm * 
m=1 


Comparing the equations (7) and (17), we find 


m (2 n) 4/2 (—1)"+1 e—naia Lea: 


A,= aria 2; aes (18) 


where 7, is given approximately by 


15 a 
ck Name Sunn ~ 
Having obtained the unknown coefficients A, substituting them into (4), 
we can find the stresses for y < 1 as 


%,=GOr+G 3) AT ("2") cos™2 , | 
pote a (19) 
ao Suan (eetantae, | 
n=l 
where 
2(_4)\n4+1 1/2 p-naja © 
as ee ie enna, Lm. (20) 


a 
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If we denote the couples acting at the ends of the shaft by T, we have 


"Oy 
T=2nG | ar=22Glyh, (21) 
0 


since y vanishes at r = 0. 

The complicated nature of the foregoing analysis is due to the abrupt 
change of the diameter along the length. If the diameter of the shaft varies 
gradually, using the expression (7), the coefficients A, may be determined 
directly so as to satisfy the boundary condition (3) by a usual approach, and 
consequently the analysis can be far more simplified. But if the change in dia- 
meter is abrupt, the solution, referred to the coordinates (7, z), often has 
singular points within the domain bounded by the cylindrical surface circum- 
scribing the shaft and the outer surface of the shaft. On the other hand, the 
series (7) must be convergent at every point within the shaft, viz. even for the 
largest value of 7. Accordingly, by means of a usual procedure, mathematical 
difficulties are encountered for the determination of A, satisfying these counter 
requirements simultaneously. 

To avoid these difficulties, the use of curvilinear coordinates is sometimes 
desirable, and the inverse method, viz. assuming first a suitable solution and 
subsequently determining the corresponding shape of a shaft, is often even 
more useful. The latter method, as is well known, considerably simplifies the 
analysis. These two methods, however, are not suitable for the cases treated in 
this paper, where the contour of the axial section of a shaft can not be repre- 
sented by a single curve. 

As an illustration of the numerical evaluations let us consider a shaft with 
grooves, the pitch and depth of which are the same as those of a Whitworth 
screw, e. g., we take 


GO13TR. BS 127. 


Assuming that C; = Cg = ... = 0, we shall find the constants Co; Gih.Ce,. Cg, 
and C, from the first five ange in (14). We thus obtain for unit 6 


Cy = 0-2952, Cy = —0-00671, C, = —0-01693, 
Cz = 0-:00301, Cy = —0-00015 . 
Substituting the numerical values for C,, into (20), we have for unit 0 
A* = —6-848 x 10-1,  A# = 7-087 x 10720," A? == 7-833,.« 105m 
—10-03.x 10-".; 


l 


Af = 8-823 x 10-“, A 
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From equations (15) and (19), we have the shearing stress on the section z = 0, as . 


= GOr+G YAEL ("2"), (r <1) 
Sins . 
=G{Or4+ 57 5? RIC, p+ Co G2 + Cz ps + Ca Qa] (22) 
-1/2 Op OP. OPs OP, \ SS 
T9=0, : 
where : 
Og, n i 
R{ 0c | EL Fane | 
‘ 0 ar’la 
Rip.) = = R[G2| = bog + 07"), 
a O”. m2 v2 72 ~ (—1)"+1 = arte 
a Op m3 7/8 ney’ = (—1)"*1 —nar'ja . 
Gal ag ra R| Acs = Canoea is a. ns © ‘ 
at y’4 ni 2 7 x* ay (—1)"t2 aria , 
Riva) = Sa gts aa ghee a a ee 


n=1 


The shearing stress on the section z = a is found from equation (19), as 


t= GOr4G SAS (1) 1, (7>=*), T,,=0. (23) 
n=l 


a 


The distributions of the shearing stress tT), on the sections, z = 0 and a, cal- 
culated from equations (22) and (23), are shown in Figure 2. It may be seen 


lal 


[unit Gy ) 
= eee 


SHresses 
SS 

= 

Bas 


S 
T 


ee Pr 
0 G20 88% Gt Se uitee 
Fig. 2 
Stresses on the sections, z = 0 and a, when a = 0-137, A= 1272 
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from the figure that in the expression for stresses (19), all terms become very 
small quantities as compared with the first term of t,,, when 7 < 0-8. Conse- 
quently, the solution for a solid shaft obtained here can be immediately 
applied to a hollow shaft, except for the case of a very thin tube whose mini- 
mum thickness is smaller than the depth of notches. The shearing stress along 
the boundary curve of the shaft can be calculated from equation (16). In our 
evaluations, however, the boundary condition (10) is not accurately satisfied, 
since only a few terms are taken in the series for y*. To estimate the error for 
the evaluations described, we shall obtain the shearing stresses, normal and 
tangential to the boundary curve of the shaft, from the equations 


i 1 J 
Th = (i + A2)22 {ht ak sii iz, al} gg (1 + A212 {[t% Ne hl) [zs Bey , (24) 
respectively, in which 


Teh GOL + haz’) + i (1 + haz')-9? RIC, g, 4+ CoG. + Coon + Cs Ga), 
Oz4s 2 


2 i) @) 0 fi 


; : ) : 0 
[egh—=—G + haz) REC, SH + iC SH iC +50, SH, 


It ay IZ 
— {1 + e-*hz cos") 
a a 


Op, 5 
R| 0c i: 3 —2 he y —nhz TZ 3 
US (eK Pt aE LCOS = 
[o,0) q — hz y 
a Oe] _ , (=D) een e N12 
R[9,], = 2 RSP] = ahe ea} n vant rot ee 
n= 
a Os] _ ** (19 2 =) cae 
BRIA - {Flee ], zl — a) + 
eS) (ny yent ga maha, naz 
jas n2 i Oe 
n=1 
g 22 hia! 
a OPs] _ TE Ar Gee 3he -) I 
Riel =< R[5'],= (12 ro eS 
lee) —nm hz 
(—1)"+1 ens NZ 
eo n® a iki a 
n=1 
4 6 h? 2/3 22 ae (ea fk we te sl : Tat 
Rip), = sq (#2 se te) tag (ME Sas) > 360 


es (—1)"+2 e—nahz nnez 


_ a + ¢os-———— 
oy, n* a 


n=1 
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It fee 
— eth?’ sin —— 
a a 


. O~ 
Rj; | = => 
0c s —27he 5 —m hz’ sbise 
Lie BNI SOE BO 
0° — hz z 
oo. AZ I —1)"+1 ens eI Tae 
Ri <2] =- . Fila ya ) ey aT 
OC Is an al" a n 
3hzz! cy (—1)Mtle-mahy | naz 
Rl pos | ke oe ca ee sin : 
Of Js a? a, n2 a 
- Os at z 2 ue 
pa) | ee zs £.ft 2 oN eee 
Rj aI, 2% (3 a2 h?z )- 25 


fore) Pe 
a —])"+1 e-nahz , NUE ‘4 Zz 
aia De kiko sin - z=1—— 
a a 
n=1 


n> a 


The shearing stresses at the boundary, calculated from equation (24), are 
shown in Figure 3, in which s is the distance from the apex along the boundary 
curve. As is shown in the figure, t,, does not entirely vanish along the curve, 
but it becomes a negligibly small quantity as compared with t,. Hence, the 
approximate solution obtained here is sufficiently accurate for practical use. 


48 


Stresses{urit Gg) 
SS S& ~ 
SS Ss 


S 
SS 


a2} 


Fig. 3 
Stresses on the bounding curve, when a = 0-187, k= 1:27. 


From (21), the couples acting on the ends of the shaft are 


T =22GC,=1:855G0. 


This shows that the stiffness of a threaded shaft is slightly larger than that of 


a shaft without grooves, the diameter of which is the same as the minimum 
diameter of the threaded shaft. 
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Zusammenfassung 


Diese Arbeit enthalt eine Methode zur approximativen Ermittlung der Span- 
nungen in einer auf Torsion beanspruchten Welle mit einem Kreisquerschnitt, 
dessen Durchmesser periodische, gegeniiber den Dimensionen der Welle kleine 
Schwankungen aufweist. Ausfiihrlich behandelt ist der Fall, wo der Durchmesser 
periodisch linear zu- und abnimmt. Fiir eine Welle, die als Modell einer Whitworth- 
Schraube aufgefasst werden kann, sind die Spannungen numerisch ausgewertet. 


(Received: December 3, 1952.) 


Zur Berechnung der Druckpunktverteilung tiber die 
Spannweite ftir Fliigel mit kleinem Seitenverhiltnis 


Von ERNst SCHULTZE, Altenrhein?) 


Bei der Berechnung der Luftkrafte an Fliigeln mit kleinem Seitenverhaltnis 
ergeben sich Schwierigkeiten, die bei Fliigeln mit grossem Seitenverhaltnis nicht 
auftreten. Fliigel mit grossem Seitenverhaltnis lassen sich mit sehr guter Nahe- 
rung durch eine tragende Wirbellinie ersetzen, deren Starke von der Mitte des 
Fliigels nach aussen stetig abnimmt. Zur Bestimmung der Wirbelstarke dieses 
gebundenen Wirbels, welche ja proportional zum Auftrieb ist, hat man die 
Bedingung, dass die von den abgehenden freien Wirbeln erzeugte Abwarts- 
geschwindigkeit gerade lings des gebundenen Wirbels mit der Abwartsge- 
schwindigkeit am Fliigel selbst ttbereinstimmt. Dies fiihrt auf die bekannte 
-Prandtlsche Integrodifferentialgleichung fiir die Auftriebsverteilung langs der 
Fliigelspannweite. Es ist aber einleuchtend, dass dabei die Bedingung der tan- 
gentialen Strémung tiber die ganze Tragflache um so weniger erfillt sein kann, 
je kleiner das Seitenverhdltnis wird (Figur 1), so dass die Auftriebsverteilung 
und der Gesamtauftrieb bei kleinem Seitenverhaltnis nicht mehr richtig heraus- 
kommen. Uber die Momentenverteilung langs der Spannweite sagt die Prandtl- 
sche Theorie nichts aus. Man muss vielmehr annehmen, dass sich der Druck- 
punkt in jedem Schnitt des Fliigels wie in ebener Strémung in 1/4 der Fliigel- 
tiefe befindet. Bei grossem Seitenverhiltnis wird diese Annahme eine gute 
Naherung geben, bei kleinem Seitenverhaltnis wird aber der Druckpunkt nicht 
mehr in 1/4 der Fliigeltiefe sein, sondern er wird sich noch langs der Spann- 


weite andern. 


1) Flug- und Fahrzeugwerke AG. 
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Im Laufe der Zeit sind einige Rechenmethoden zur genaueren Berechnung 
von Auftriebsverteilung und Druckpunktverteilung aufgestellt worden. Eine 
bekannte Methode ist diejenige von V. M. FaALKNER}). Er ersetzt den Fligel 
durch ein Gitter von Wirbellinien, das 84 bzw. 126 Gitterfenster enthalt. Damit 
bekommt dann die Auftriebsverteilung iiber die Spannweite einen treppen- 
férmigen Verlauf. Obwohl diese Methode einen gewissen Fortschritt gegentiber 
PRANDTL brachte — sie lasst sich auf beliebig kleine Seitenverhaltnisse und 


v gebundener Vv 
_—  Wirbel 
ae 
free 
Wirbel 
Grosses Seitenverhaltnis Kleines Seitenverhaltnis 


Fig. 1 
Zur Erklarung des mit abnehmendem Seitenverhaltnis wachsenden Einflusses der durch die freien 
Wirbel induzierten Kriimmung der Stromlinien. 


beliebig grosse Pfeilung anwenden —, so bedeutet der Ubergang von der ste- 
tigen Auftriebsverteilung der Prandtlschen Theorie zu einer unstetigen Auf- 
triebsverteilung einen Riickschritt. Deshalb sind andere Methoden entwickelt 
worden, die als direkte Verallgemeinerungen der Prandtlschen Traglinien- 
theorie bezeichnet werden kénnen. 

Einen ersten Schritt in dieser Richtung macht die Theorie von WEISSINGER®). 
Er legt die tragende Linie in 1/4 der Fliigeltiefe und erfiillt die Bedingung der 
tangentialen Strémung in 3/4 der Fliigeltiefe. Die dadurch entstehende Integro- 
differentialgleichung fiir die Zirkulationsverteilung ist nur wenig komplizierter 
als die Prandtlsche und kann ebenfalls nach MuLtuoppP gelést werden. Die sich 
daraus ergebenden Auftriebsverteilungen stimmen auch bei kleinen Seiten- 
verhaltnissen und grosser Pfeilung in den meisten Fallen sehr gut mit dem 
Experiment iiberein (nach NACA Report Nr. 921, 1948). Auch die Differenz 


1) V. M. FaLkner, British Reports and Memoranda, Nr. 2591. 
*) J. WeIssincErR, Nat. Advis. Comm. Aeron. TM. 1120 (1947). 
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gegentiber der nach FALKNER berechneten Auftriebsverteilung ist gering, so 
dass es sich nicht lohnt, zur Berechnung der Auftriebsverteilung die miihsamere 
Falknersche Methode zu beniitzen. Hingegen sagt auch die Weissingersche Theo- 
rie tiber die Druckpunktlage nicht mehr aus, als dass der Druckpunkt eben dort 
liegt, wo sich der tragende Wirbel befindet, naémlich in 1/4 der Fliigeltiefe. 

Um eine Momentenverteilung zu bekommen, ist es deshalb nétig, mehr als 
einen tragenden Wirbel anzunehmen. So hat N. ScHo1z!) vier solche Wirbel- 
linien angenommen. Er teilt den Fliigel der Lange nach in vier gleich breite 
Streifen ein. Die Wirbellinien befinden sich in 1/4 der Tiefe jedes Streifens und 
die Linien, langs denen die Abwindbedingung erfiillt wird, in 3/4 der Tiefe 
jedes Streifens. ScHoiz berechnet damit die Auftriebs- und Druckpunktver- 
teilungen von rechteckigen Tragfliigeln mit den Seitenverhaltnissen 10, 6, 3 
und 1. Es ist natiirlich klar, dass man um so genauere Resultate bekommt, je 
mehr Wirbellinien man annimmt. Anderseits wird man im Interesse des klein- 
sten Arbeitsaufwandes ein Minimum an Wirbellinien anstreben. 

Im folgenden wird nahegelegt, wie man mit zwei Wirbellinien und zwei 
Linien der tangentialen Strémung, die sich je in geeigneter Fliigeltiefe befinden, 
ebenso gute Ergebnisse erhalten kann wie mit der obigen Anordnung fiir vier 
Wirbellinien. 

Der Grund, warum die Weissingersche Theorie so gute Auftriebsvertei- 
lungen liefert, kann an Hand der ebenen Strémung um diinne Profile erklart 
werden. Es liegt die Vorstellung zugrunde, dass bei einer diinnen ebenen Platte 
von endlichem Seitenverhaltnis die abgehenden Wirbel eine Kriimmung der 
Str6mung erzeugen, so dass sich kein Fliigelschnitt mehr in ebener Strémung 
befindet, dass dies aber auf dasselbe herauskommt, wiewenn jeder Fliigelschnitt, 
nachdem er entsprechend verwélbt worden ist, sich als in ebener Stro6mung 
befindlich betrachtet werden kann. Da man die ebene Strémung um gewdlbte 
Profile gut beherrscht, bleibt also nur noch das Problem der Bestimmung 
dieser fiktiven Verwolbung, die durch die abgehenden Wirbel induziert wird. 
Wir gehen von der Birnbaumschen Theorie der ebenen Strémung aus, in wel- 
cher drei Normalfalle betrachtet werden (Figur 2). Diese entsprechen den ersten 
drei Gliedern der Entwicklung der Wirbelbelegung in eine Fourier-Reihe. Die 
iibrigen Glieder dieser Reihe ergeben keinen Beitrag zum Auftrieb und zum 
Moment. Ersetzt man die stetige Wirbelbelegung eines Profils in ebener Str6- 
mung durch einen einzigen Wirbel von der Starke I’, so muss, damit er den 
richtigen Auftrieb liefert, die Bedingung 


T= tq (1) 


erfiillt sein. Der Abstand dieses Einzelwirbels von der Vorderkante sei ¥,, und 
der Abstand der Stelle, wo die tangentiale Stromungsbedingung erfiillt werden 


1) N. Scuorz, Ing.-Arch. 18, 2 (1950). 
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Ebene Platte Kreisbogenprofil S-Schlagprofil 


0 Cc Ae 
5 M4 
Cg 2s Ogi mee ‘ 
x 
w=vae 42-29 
c[2 : 
Die Momente sind auf die Vorderkante (x = 0) bezogen und positiv, wenn sie schwanzlastig 


wirken. Der Abwind w ist nach unten positiv. v = Fluggeschwindigkeit, c = Fligeltiefe. 


Fig. 2 


Die drei Birnbaumschen Normalfalle eines ditinnen Profils in ebener Str6mung. 


soll, von der Vorderkante sei x». Die beiden Unbekannten x, und x, sind nun 
so zu bestimmen, dass sowohl die ebene Platte fiir sich als auch das Kreis- 
bogenprofil fiir sich in Hinsicht auf den Auftrieb durch den Einzelwirbel 
ersetzt werden kénnen. 
Ebene Platte: 

Auftriebsbedingung: 


EEE 


Abwindbedingung: . 


daraus durch Division 


%y—M= >. (2) 
Kreisbogenprofil: 
Auftriebsbedingung : 
l= a vcB, 
Abwindbedingung: 
is ap) ps x . 
Fae = 2% w(%q) = 2x vB (a5 —1); 
daraus durch Division 
(% — %) (45-1) =4 (3) 
Die Auflésung der beiden Gleichungen (2) und (3) ergibt 
3 1 
XO geo. er ae (4) 
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Daraus folgt 

Satz 7. Ein diinnes Profil, das aus einem Kreisbogen beliebiger Wolbung 
bei beliebigem Anstellwinkel besteht, kann hinsichtlich des Auftriebes durch 
einen in 1/4 der Fliigeltiefe liegenden Einzelwirbel ersetzt werden, wenn die 
tangentiale Strémungsbedingung in 3/4 der Fliigeltiefe erfiillt wird. 

Der grosse Erfolg der Weissingerschen Theorie der Auftriebsberechnung 
liegt also in der Beriicksichtigung der durch die abgehenden Wirbel induzierten 
Kriimmung der Stromlinien, wahrenddem die Prandtlsche Theorie nur den 
durch die abgehenden Wirbel induzierten Anstellwinkel beriicksichtigt. Hin- 
gegen versagt die Weissingersche Theorie bei der Berechnung des Momentes, 
da der Druckpunkt einer einzelnen Wirbellinie immer an der Wirbellinie selbst 
liegt, in diesem Fall in 1/4 der Fliigeltiefe. Fiir das Moment wird also nicht die 
induzierte Kriimmung der Stromlinien, sondern nur der induzierte Anstell- 
winkel beriicksichtigt. 

Der Weg zur Berechnung der Momentenverteilung ist nun aber vorge- 
zeichnet: Wir miissen zwei Wirbellinien nehmen. Ihre Wirbelstarken seien mit 
I’, und I’, bezeichnet und ihre Abstaénde von der Vorderkante mit x, und 1». 
Die Auftriebsbedingung liefert 


P+ T,= [ye dx= 7 cg (5) 
0 
und die Momentenbedingung, welche nun auch erfiillt werden soll, 
; Uc? 
mT, + %T,= | yo) xdx=—* en. (6) 
0 


Die Abwindbedingung soll wiederum im Abstand x, hinter der Vorderkante 


erfillt werden: 
rT; I, 
A —— Leal Xy — Xe 


= 27 W(X) - ' (7) 


Im Falle der ebenen Platte liefern diese drei Gleichungen: 
Sa pe ks =Tvcta, (8) 


ie Sa ee 
Xo — #y Xo — Xe 


= ES C 


Durch Elimination von I’, und I, gelangt man zu der folgenden Gleichung, 
in welcher zur Abktirzung die Grdssen &) = xo/c, &, = x4/¢, & = %/¢ eingefiihrt 


wurden: 


54+ £-&- &=2 (6 - &) £0 - (9) 
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Dies ist eine Gleichung fiir die drei Unbekannten &, —, und é,. Eine weitere 


Gleichung erhalt man aus dem Fall des Kreisbogenprofils, welcher zunachst 
aus (5), (6) und (7) liefert 
ee or ete 7 + VB, 
ty ly + tele = P, (10) 
ra ne gos = 2nv8 (35-1) 
Nach Elimination von J’, und J: 
5 +8 —&:- & = 4 (20-1) Eo— A) Ge &)- (11) 


Eine dritte Gleichung zwischen den drei Unbekannten & , &, und &, ergibt sich — 


aus dem Fall des S-Schlagprofils 


r,t Se” SOR | 
aig 
m0, + mI, =-y ved, (12) 
T; I eee % \2 Xo | 
Hy — Hy rg Ho—%y 2 ge [-2 Ga = o/2 -1] E : 
durch Elimination von I’ Sata 
1 = 16 (1 — 85) + 8 5) (So — 4) (Fo — 5) - (13) 


Die Auflésung der drei Gleichungen (9), (11) und (13) nach den drei Unbekann- | 


ten geschieht folgendermassen: Zuniachst gilt: 
(Fo — &1) (0 — &) =o — még tm, 


wobei 


ra 


n=&,+ und m=6,é, (14) 
gesetzt wurde. Sodann dividiert man (9) durch (11). Dies gibt 
8&@—4& (1+2n)—2+6n=0. (15) 
Subtrahiert man hievon das Vierfache der Gleichung (9), so bleibt: 
Sm—2n+1=0. (16) 
Die Division (9) : (13) ergibt 


64 & — 10 & (3+4n) —8& (1—8n)+1-—8n=0. 
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Subtrahiert man hievon das (8 &, — 2)-fache der Gleichung (15), so folgt 


4% = 3, 
also zusammen mit (16) 
, Psy ees 
4? 16° 


Wegen (14) sind &, und &, die beiden Lésungen der quadratischen Gleichung 


3 1 
16 0 ; 
also 

5 


ae 
BV Bie 9 0955! le = 0 I 


5 - = 0,6545. (17) 


Aus Gleichung (15) ergibt sich durch Einsetzen von m auch eine quadratische 
Gleichung fiir &,: 


48-5842 =0 


mit den beiden Lésungen: 


(18) 


Daraus folgt 

Satz 2. Ein diinnes Profil in zweidimensionaler Str6mung, das aus den An- 
teilen: ebene Platte, Kreisbogenprofil, S-Schlagprofil besteht, kann hinsichtlich 
des Auftriebes und des Momentes in exakter Weise durch zwei gebundene 
Einzelwirbel ersetzt werden, deren Abstande x, =cé, und x, =cé&, von der 
Vorderkante durch die Formeln (17) gegeben sind, wenn die tangentiale Str6- 
_mungsbedingung in den zwei durch die Formel (18) gegebenen Linien im Ab- 

stand x) = c €) von der Vorderkante erfiillt wird. 

Wie sich auf dem Satz 1 die Weissingersche Theorie der Auftriebsberechnung 
aufbaut, so liefert der Satz 2 eine Methode zur Berechnung der Momentenver- 
teilung. In dieser Methode werden also zwei gebundene Einzelwirbel angenom- 
men, die in den durch die Formeln (17) gegebenen Fliigeltiefen liegen. Die von 
diesen gebundenen Wirbeln und den abgehenden freien Wirbeln induzierte 
Abwartsgeschwindigkeit wird dann lings zwei Linien, die in den durch die 
Formel (18) gegebenen Fliigeltiefen verlaufen, gleich der durch den Fligel 
gegebenen Abwartsgeschwindigkeit gesetzt. Dies ergibt ein System von zwel 
Integrodifferentialgleichungen fiir die beiden unbekannten Zirkulationsver- 
teilungen, welches sich auf dieselbe Weise durch ein lineares algebraisches 
Gleichungssystem ersetzen lasst wie die Weissingersche Integrodifferential- 
gleichung, nur mit dem Unterschied, dass die Anzah] der Unbekannten jetzt 


o>? 
doppelt so gross ist. Auf diese Weise wird dann also nicht nur der durch die 
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rdumliche Strémung induzierte Anstellwinkel, sondern auch die induzierte 
Kreisbogen- und S-Schlagkriimmung der Stromlinien berticksichtigt. Auch 
Tragflachen, die selbst Kreisbogen- und S-Schlagkriimmung aufweisen, lassen 
sich behandeln, wenn man nur an den Stellen, wo die tangentiale Strémungs- 
bedingung erfiillt wird, die entsprechende Profilneigung einsetzt. 


Summary 


After having shown that the PRANDTL and WEISSINGER theories for the cal- 
culation of the lift distribution of wings of finite span fail to give the spanwise 
distribution of the center of pressure, a new model replacing the wing, consisting 
of two bound vortices, is proposed for the calculation of the lift- and moment 
distribution. The positions of the two bound vortices and of the two lines for 
satisfying the downwash condition are determined—in a straight-forward 
manner as an extension of the WEISSINGER model—so that the first three terms 
of the FourtER development of the chordwise vorticity distribution give exact 
values of lift and moment at every spanwise station. 


(Eingegangen: 20. Dezember 1952.) 
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Rechnen mit komplexen Zahlen auf einer mechanischen 
Multipliziermaschine 


Von WALTER Bo6scu, Ziirich?) 


Die numerische Berechnung des resultierenden komplexen Widerstandes von 
zwei parallelgeschalteten komplexen ,Widerstanden auf einer gew6hnlichen mecha- 
nischen Multipliziermaschine kann dann erfolgen, wenn der Ausdruck: 


Se ee 


(wobei 3, = Fh, + 7 X, und 3,= R, + 7X, ist) iibergefiihrt wird in die Gleichung: 
31 (RE + X$) + Be (RT + X73) 


Sp TR + Ra) + +X)" 

Gegeniiber der bis heute allgemein iiblichen Berechnungsmethode mit Leit- 
werten und trigonometrischen Funktionen ergibt sich in der Praxis eine héhere 
Genauigkeit und eine grosse Arbeitszeitreduktion, speziell bei Netzberechnungen. 


(Eingegangen: 14. Februar 1953.) 


1) Postfach Ziirich 23. 


Vol. IV, 1953 Kurze Mitteilungen — Brief Reports — Communications bréves Zi 


Kleine Schwingungen dynamischer Systeme 


Von Paut Ernés, Ziirich!) 


Einleitung 


Die Bewegung eines holonomen konservativen dynamischen Systems mit 7 
Freiheitsgraden um seine Gleichgewichtslage oder um seine gleichférmige Bewe- 
gung wird durch die den Anfangsbedingungen gentigenden Losungen der Hamil- 
tonschen Differentialgleichungen 


dn =- on q yee 
rer gs 04K 


bestimmt, wobei ¢;, bzw. ¢,,.; die verallgemeinerten Lage- bzw. Impulskoordinaten 
H die Hamiltonsche Funktion des Systems bezeichnen. 

Beschrankt man sich auf kleine Bewegungen des Systems, so kann die 
Hamiltonsche Funktion als quadratische Form der Variablen qj, ..., dz, ange- 
geben werden. Dabei kénnen gyroskopische Terme der Form dp Gn4; (k, 2 < ”) 
auftreten, die es verunméglichen, die Stabilitat durch eine simultane Haupt- 
achsentransformation der kinetischen und potentiellen Energie zu untersuchen. 

Es gelang Weierstrass (1]?), unter der Annahme, dass Sf positiv definit ist, 
die Stabilitat der Gleichgewichtslage bzw. gleichfoérmiger Bewegung aufzuzeigen, 
und die explizite Lésung der Hamiltonschen Differentialgleichungen in Form einer 
Linearkombination harmonischer Schwingungen anzugeben. 

Die WeierstraBsche Methode wurde auch von WHITTAKER 2] tibernommen. 
Die neuere Arbeit von FrAzER, DUNCAN und CoLLar [8] verweist auf WHITTAKER. 
Besonders die Uberwindung der beim Auftreten mehrfacher Wurzeln der charak- 
teristischen Gleichung von & entstehenden Schwierigkeiten (die LAGRANGE [6], 
LapLace({7] et al. zu Fehlschliissen verfiihrten) gibt zu einer langwierigen mathe- 
matischen Behandlung Anlass. 

Im ersten Abschnitt soll gezeigt werden, wie sich durch Anwendung von Me- 
thoden der linearen Algebra eine kurze Behandlungsweise ergibt, die sich von der 
WeierstraBschen unterscheidet. 

Im zweiten Abschnitt werden diese Methoden der linearen Algebra auf den 
WeierstraBschen Lésungsgang iibertragen, und der Zusammenhang mit dem im 
ersten Abschnitt beschriebenen Vorgehen hergestellt. 

Im dritten Abschnitt wird ein Hilfssatz der Matrizentheorie bewiesen, der 1m 
ersten Abschnitt Verwendung fand. 


fc ee) 


Bezeichnungen 


Falls nichts anderes vermerkt wird, bezeichnen kleine fette Buchstaben 
(2 n)-dimensionale Vektoren (Kolonnenmatrizen). So sei zum Beispiel g der Vektor 
mit den Komponenten 4,(é), .--, Zen(t)- Grosse fette Buchstaben bezeichnen qua- 


dratische Matrizen vom Grade 2, zum Beispiel 
OE 
@=(_ Fo): 


1) Physikalisches Institut der ETH. 


2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 219. 
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wobei E, O die Einheits- bzw. Nullmatrix vom Grade nu bedeuten. Es ist 
Q* = Q:''= —-Q. Der * bezeichnet die transponierte Grésse. Die Gesamtenergie 
des Systems schreibt sich 


1 F 
H=—>+q@q*Hq mt H*=H. 
Ferner soll eine charakteristische Matrix zE — A als eine solche mit lauter 
linearen Elementarteilern bezeichnet werden, wenn 2; eine «,-fache Wurzel der 
Ls 
charakteristischen Gleichung f(z) = |zE— A|=Oist, und 7 (z— z,)**~1 als grésster 
k=1 


gemeinsamer Teiler aller (2 » — 1)-reihigen Unterdeterminanten derselben auftritt. 
Das Produkt erstreckt sich dabei nur iiber die y voneinander verschiedenen Wurzeln. 


tb 
Fasst man die 2 2 Hamiltonschen Gleichungen zu 
g-QHq (1.1) 
zusammen, so kann g mittels der zeitabhangigen Matrix B(f) als 


oo 
(¢— 4)” 
q = BY) g@= DB, a7 
=0 
angesetzt werden, mit g = g(¢)) den Anfangsbedingungen Rechnung tragend. 
Um die Beiwerte B,, der nach der Zeit in Potenzreihe entwickelten B(é) zu be- 
stimmen, wird g und 


$= BY) ¢ =D) Bae d (1.3) 
in (1.1) eingesetzt: se 
booed m oS m 
Bren ("4 on Sa, 8 
Der Koeffizientenvergleich ergibt 
und damit sigs i: — 
.y m (E— to)" 6 i Qui-t,) 2 
q 2, (2 FS nn q- (1.5) 


Um diese Lésung in der iiblichen Form darzustellen, verwendet man die 
Eigenschaft von Q H, dass ihre charakteristische Matrix z E — Q A lauter lineare 
Elementarteiler und paarweise konjugiert imaginare charakteristische Wurzeln 
2; = +1 (x; reell, positiv) besitzt, da der grésste gemeinsame Teiler aller ihrer 
(2 — 1)-reihigen Minoren mit dem der aquivalenten Matrix 


H>? H(z E— QH) Hl =2E-— H'2Q HP (1.6) 


1) Dies ist eine unmittelbare Folge davon, dass die Determinante eines Produktes von Matrizen 
das Produkt der Determinanten der Faktoren ist. Siehe zum Beispiel [4]. 
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(H*!" Q Ht!*)* = —H)? Q HH” somit H1/? Q H!!? schiefsymmetrisch. H1!? Q H12 
hat die obenerwahnte Eigenschaft [3], folglich auch QH. 

Hilfssatz. Hat die charakteristische Matrix einer Matrix P lauter lineare 
Elementarteiler und ist g(z) eine regulare eindeutige analytische Funktion, so gilt 


Z— By 
sim = Sten (Fe)... a 
wobei 2;, ..., 2, die voneinander verschiedenen Wurzeln von f(z) = |z E—P|=0 


sind. 
Wendet man dies auf Q Han, so wird 


4}2 . . 
mo ao Ut, ee Zt; 
Gi eixglt o( she) ) | g—ixg(t—to) (, Cj ) 3 1. 
S ZE— OH )s=ix; A Te): Elo! sae aaa 


Bezeichnet man 


Z— 1%; _/( @—t2) @E=Q HH) _ ; 
Cerner (2) [eee Py) i) 
(F;, Fj reell), so ist 
7/2 
q=2 be cosx,; (f — to) — Fj sinx, (¢ — | g. (1.10) 
7a 
Ide 
g—QHq (et) 
wird auch durch den Ansatz 
— b Zz e2(t—to) dz Za 
q (2.1) 
6; 


gelost, sofern als Integrationsweg C ein Kreis gewahlt wird, der samtliche Wur- 
zeln von f(z) = |zE—QH|= 0 enthalt. Setzt man diesen Ansatz fiir q und 
dessen Ableitung nach der Zeit fiir g in (1.1) ein, so kommt man auf die Bedingung 


f (zg E— QA) b(z) e?"-%) dz=0, (2.2) 


die nur erfiillt sein kann, wenn (z E— QH) &(z) eine in C regular-analytische 
Vektorfunktion 


[oe) 
(2 E- QH) b(2) = 3%, = 9) (2.3) 
k=0 
ist. 
In einer ringférmigen Umgebung des Integrationsweges lasst sich (2.3) nach 
6(z) aufldsen und liefert 


b(z) = (*#E— QA)" gz). (2.4) 
Somit ist die Losung: ta 
4 ge a(t to) 2.5 
q f BE = OH e dé. (AR) 
C 


1) Die Annahme [1], [2], dass mit Ausnahme von gp alle Entwicklungskoeffizienten Null sind, 
ist nicht erforderlich. 
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Beachtet man die im ersten Abschnitt hergeleitete Tatsache, dass ( E— Q H)~* 
nur einfache, konjugiert imaginare Pole hat, so ergibt die Anwendung des Resi- 
duensatzes auf dieses Integral den Ausdruck (1.10). Die Beriicksichtigung com 
Anfangsbedingungen fihrt auf 
2%ty(z) = q. 
Stellt man sich die Frage, welcher logische Weg zum Ansatz (2.1) bzw. (2.5) 
fiihrt, so ergibt sich die Antwort aus dem Vergleich von (1.5) und (2.5): 


yi e@2(t—t) 
H(t=t,). = - at ae 2s 
a rat PEE= OHS “= 


In der Tat kann eine eindeutige analytische Funktion einer Matrix A in all- 
gemeinster Weise durch das Cauchysche Integral 


1 g(z) dz 
2uat J zE—A 
Cc 


cya) = (2.7) 


dargestellt werden [5]. Es liegt daher nahe, den Lésungsvektor als ein solches — 
Cauchy-Integral anzusetzen, wobei die Zeit als Parameter auftritt. 

Die obigen Betrachtungen liefern auch ein Beispiel dafiir, dass die spektrale 
Zerlegung einer Matrix vollkommen gleichbedeutend ist mit der Anwendung des © 
Residuensatzes auf ihr Cauchysches Integral. 


3. 


Beweis des Hilfssatzes. Ist P eine Matrix vom Grade », deren charakteri- 


stische Matrix lauter lineare Elementarteiler besitzt, und 
P= pa a; P; (3.1) 
t 
ihre spektrale Darstellung, mit i 
E=)'P;, (3:2) ; 
t ’ 
P,; Pj = 055 P; (3.3) 
} 


gehorenden Eigenraum langs des komplementaéren Raumes projizieren, dann 


wobei die P; Operatoren darstellen, die jeden Vektor auf den zum Eigenwert 2; 
besitzt (¢ E — P)~* die spektrale Darstellung ; 


Pa (@QE=— Peo 1 
zE—P fe ae ay Pee (3.4) 
und es ist 
: Z— & j 
7 ey el (35) 


Ist z, eine «,-fache Nullstelle von f(z), so ist sie — geméass Voraussetzung — eine 
(%;, — 1)-fache jedes Elementes von (z E— P)*%, da die letzteren gerade die 
(m — 1)-reihigen Minoren von z E — P sind. Folglich existiert obiger Limes und 
ist von Null verschieden. (Fiir Matrizen mit nichtlinearen Elementarteilern exi- 


stiert dieser Limes nicht und somit auch keine spektrale Darstellung im obigen 
Sinne.) 
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Ist g(z) eine regulare eindeutige analytische Funktion, so ldsst sie sich in eine 
Laurent-Reihe entwickeln, und es gilt wegen (3.13): 


g(P) = 2 8 (2;) P; - (3.6) 


Zum Schlusse danke ich den Herren Prof. Dr. E. StreFeL*und PD. Dr. 
H. RUTISHAUSER herzlich fiir ihr f6rderndes Interesse an der Arbeit. 

Mein Freund Hans WIDMER hat mir in verdankenswerter Weise bei der Text- 
korrektur beigestanden. 
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Summary 


The Hamiltonian differential equations for small displacements of a dynami- 
cal system with a finite number of degrees of freedom are solved by use of the 
methods of linear algebra. In this way the stability of the equilibrium configura- 
tion or uniform motion in the case of positive definite total energy is demonstrated. 
Especially gyroscopic terms may appear in the hamiltonian, making impossible 
the usual proof of stability by means of transformation to normal coordinates. 
This method is then applied to the Weierstrass procedure of solution. 

For use in the above intefration, the spectral representation of a matrix with 
linear elementary divisors is given through the adjoint of its characteristic matrix. 


(Eingegangen: 16. Oktober 1952.) 


Uber eine Gefahr, die durch Ruderausschlage bei schallnahen 
‘Geschwindigkeiten entstehen kann 


Von EuGeN Gruscuwitz, Altenrhein?) 
Es ist bekannt, dass die Ruderwirkung bei Flugzeugen mit zunehmender 


Machscher Zahl M zunichst ansteigt, so wie es die Prandtlsche Regel veriangt. 
Von einem gewissen Wert von M an, der natiirlich von verschiedenen Einfliissen 


1) Flug- und Fahrzeugwerke AG. 
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abhangt, kann sie indessen mit weiter zunehmendem MW sehr rasch absinken und 
Werte erreichen, die nur noch einen Bruchteil derjenigen bei kleinen Werten von 
M darstellen. Zum Beispiel hat sich aus einer Messung von O. KNAPPE?) an einem 
Fliigel mit 30% tiefem Ruder in ebener Str6mung der in Figur 1 dargestellte 
Verlauf der Ruderwirkung, ausgedriickt durch den Wert 0c,/0f, iiber M ergeben. 
Dieser Zusammenbruch der Ruderwirkung ist auf Stromungsablésung am Ruder, 
die durch Verdichtungsst6sse verursacht wird, zuriickzufiihren. 


rl T 


|h cz 


Qa 05 i 
Fig. 1 
Oc,/0B als Funktion der Machschen Zahl. tp/t = 0,3. 
(Nach KNApPE.) 


Um bei Machschen Zahlen nahe 1 eine betrachtliche Ruderwirkung zu er- 
reichen, benétigt man also grosse Ruderausschlage. Mit den bei den meisten 
Hochgeschwindigkeitsflugzeugen vorhandenen Kraftsteuerungen sind solche Aus- 
schlage trotz der grossen Rudermomente auch erreichbar. Erfolgt nun eine ziem- 
lich rasche Verminderung der Fluggeschwindigkeit, wie sie etwa beim Abfangen 
aus dem Sturzflug, insbesondere beim Gebrauch von Sturzflugbremsen oder 
auch beim Ubergang aus dem Horizontalflug in einen steilen Steigflug zu erwar- 
ten ist, und wird wahrend dieses U berganges der grosse Ruderausschlag vermége 
der Kraftsteuerung beibehalten, so kann die Ruderwirkung infolge der Abnahme 
von M sehr rasch auf ein Mehrfaches ihres anfanglichen Wertes steigen. Die zum 
Abfangen nétige Héhenleitwerkslast liegt im allgemeinen nicht sehr weit unter 
der sicheren Last, fiir die die Festigkeit des Leitwerkes dimensioniert ist. Durch 
den eben beschriebenen Effekt kann die Héhenleitwerkslast auf ein Mehrfaches 
der Abfanglast steigen. Es besteht also die Gefahr eines Leitwerksbruches oder, 
infolge zu heftigen Abfangens, die Gefahr einer Uberschreitung des zulassigen 
Lastvielfachen und damit die eines Fligelbruches. Méglicherweise lassen sich 
einige bisher ungeklarte Unfalle so verstehen. 

Es ist denkbar, dass ein ahnlicher Effekt auch ohne Anderung der Geschwin- 
digkeit entstehen kann. Die Lage des Abfalles der Ruderwirkung iiber M hangt 
vom Anstellwinkel des Leitwerkes ab. Eine Anderung dieses Winkels kann daher 
die gleiche Wirkung haben wie eine Anderung von M. Wenn zum Beispiel ein 
Ruderausschlag beim Anstellwinkel a, am Punkte 4 der Figur 1 erfolgt und der 


Anstellwinkel sich von «, auf «, andert, so springt die Ruderwirkung auf den 
Wert beim Punkt B. 


1) O. Knappe, Schnellkanalversuche an einem symmetrischen Klappenjliigel, Jb. dtsch. Luft- 
fahrtforsch. 1, 96 (1941). 
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Auf alle Falle stellen grosse Ruderausschlage bei hohen Machschen Zahlen, 
zumindest am Héhenleitwerk, eine Gefahr dar. Str6mungen mit Verdichtungs- 
stossen und Grenzschichtablosung sind der Berechnung nicht zuganglich und 
sind durch sehr kleine Ursachen beeinflussbar. Das Absinken der Ruderwirkung 
kann bald mehr und bald weniger stark und auch mehr oder weniger steil se: 
und man kann sich im Einzelfall auf Messungen nicht verlassen. Man vermeidet 
diese Gefahr am einfachsten, indem man oberhalb einer gewissen Fluggeschwin- 
digkeit an Stelle grosser Ruderausschlage entsprechende Flossentrimmungen ver- 
wendet. 

Summary 


It is called attention to the danger for airplanes resulting from the behaviour 
of the rudder efficiency at overcritical Mach numbers. Windtunnel tests show 
that in many cases at transsonic Mach numbers a sharp decrease of the rudder 
efficiency is possible. Author shows that as a consequence of this behaviour 
under certain conditions forces may be produced which are ample to destroy the 
aircraft. 


(Eingegangen: 16. Februar 1953.) 


Growth of Boundary Layer on a Rotating Sphere 


By Swami Dayar NIGAM 
and KUMANDUR SRINIVASA IYENGAR RANGASAMI, Kharagpur?) 


1. Introduction 


HowarTH?) has investigated the problem of the flow engendered by a sphere 
rotating uniformly about a diameter in otherwise undisturbed fluid. From the 
physical considerations we should expect an inflow at the poles and an outflow 
in the equatorial plane. The flow should be symmetrical about the equatorial plane, 
and the sphere may be considered to be made up of two hemispheres joined 
smoothly at the equator. HowarrH finds that the boundary layers originate at 
the poles on the two hemispheres and develop towards the equator where they 
impinge on each other. Howartu’s solutions are not valid near the equator. 
They do not give any information about the flow in the equatorial plane, He 
concludes on the basis of his solutions that inflow will occur over a large part 
of the surface, the outflow necessary to maintain continuity will be confined to 
the vicinity of the equatorial plane. His solutions give an inflow at the poles, but 
they do not give any outflow in the equatorial plane because they cease to be 
valid near the equator. 

In the present note we have discussed the growth of motion, in the earlier 
stages of its development, caused by a sphere which at the time ¢ = 0 is suddenly 
made to rotate with a constant angular spin Q about a diameter in fluid otherwise 
undisturbed. The solutions have a serious limitation in that they give initial 
motion only. They give no information regarding the time after which the steady 


state is established. 


1) Indian Institute of Technology, Kharagpur, India. 
2) L. Howarth, Phil. Mag. 42, 1308-1315 (1954). 
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2. Equations of motion 


We shall use spherical polar coordinates 7, 6, g with y measured radially out- 
wards from the centre of the sphere, 6 measured from the axis of rotation, and p 
the azimuth. If w, u, v represent the components of velocity in the directions 
y, 0, m increasing the boundary layer equations can be deduced from HOWARTH?) 
by placing K, = 0, kK, (1/a) ctg@, where a is the radius of the sphere. They are: 


Ou . lu ou Ou v2 O7u 


-——-—. Se ——_— — —— —— ¥ 1 
Ot t a& 00 hes Or a asctg ” ort’ (1) 
dv u Ov Ov uv 07u 
Lieeoter 2 bp Se = ae! 2 
FY ah or Pa ee ” Ore? (2) 
OR fig? Ow uU 
j ee ees = 0%. 3 
la 00 aS or : a vee (3) 
We now set 
u = 2? asin6@ cos@ t f(n) , (4) 
v= QDasind g(n) , (5) 
w= —2? 2 y1/2 48/2 (3 cos? 6 — 1) h(n) , (6) 
where 
¥—a 
12 (vy d)ik (7) 


During the early stages of motion when ¢ is small (or in boundary layer theory 
terminology: when the thickness of the boundary layer is small), we may neglect 
the terms in the equations of motion containing higher powers of ¢. Therefore, 
omitting terms of order ¢? in the equations of motion, we get to a first order of 
approximation the following equations: 


f+ 2nf —4f=—4 8%, (8) 
g’+ 2g’ =0, (9) 
f=h’, (10) 


where a dash denotes differentiation with respect to 7. 
These equations are the same as the equations in the problem of boundary 
layer growth over an infinite rotating disk (N1GAM)?). 


3. Solutions of the equations 


The solutions of these equations satisfying the boundary conditions u = 0, 
v= Qasin6, w=0 on the sphere y =a; and u=v=0 at ¥=&, are 


g = [1 —erfy] = eric n j (11) 
2 
f= = [(1 + 2?) erfcn — 2 2-1/2 » e-"")] — 2 (m-1/2 e-1" — 9 erfc “7 oe (12) 
2 8 2 w—-1/2 2 
b= (3.9 + 2m) erfem — 2 2/8 (1 + 9%) em) — 2 9 (a4? e- — n erfen)? 
Vids a 22 = Bia Car 
i FSi} e—™ erfc n +- 3 ly erfc (/2 n) + aur (= = att 1) (13) 


L. Howartn, Phil. Mag. 42, 289-248 (1951). 


2) 
*) S.D. Nigam, Quart. Amer. Math. 9, 89-91 (1951). 
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u vanishes both on the axis of rotation and in the equatorial plane; v vanishes 
on the axis of rotation; w is negative so long as 3cos*@—1> 0, it vanishes for 
S1c0s20 = 11-6, 6 = 54° 45" in the upper hemisphere and §6=125° 15’ in_the 
lower hemisphere, and is positive when 3 cos?@ —1< 0. This shows that there 
is a radial inflow towards the sphere for 0 < 0< 54° 45’; this radial flow dis- 
appears on the cone @ = cos~} (1//3) and changes into a radial outflow for 
54° 45’ << § < a/2. Hence there is an inflow at the poles and an outflow near the 
equator. 

It is highly probable that the cone 0 =cos~1(1///3) may viden towards the 
equator but these solutions do not give any such information. 

A stream function may be defined by the equations 


Uu i oO y Box i Ow 


a sinO or 


‘sind 00 
whence the stream function may be expressed 


y = 2 Q2 13/2 y1!2 cos6 sin? 6 h(y) . 


The flow functions. 


Zusammenfassung 


Die Gleichungen der Grenzschichttheorie wurden angewendet, um das An- 
wachsen der Grenzschicht an einer in ruhendem Wasser rotierenden Kugel zu 
studieren. Die Lésungen sind nur fiir kleine Zeitintervalle giiltig. Die Grenzschicht 
entsteht an den Polen und verstarkt sich nach dem Aquator zu. Im Gegensatz 
zu der Losung von HowartH fanden wir eine Einstroémung an den Polen und 


einen Abfluss am Aquator. 


(Received: December 29, 1952.) 
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Varia — Miscellaneous — Divers 


International Congress of Mathematicians 1954 


The International Congress of Mathematicians 1954 will be held in Amsterdam 
from September 2nd to September 9th under the auspices of ‘‘ Het Wiskundig 
Genootschap”’ (The Mathematical Society of the Netherlands). The Organizing 
Committee has invited a number of outstanding mathematicians to deliver 
one-hour addresses, hoping that in this way a survey of the recent development 
in the whole field of mathematics may be furnished. There will be seven sections, 
viz :— 

(1) Algebra and Theory of Numbers; 

(2) Analysis; 

(3) Geometry and Topology; 

(4) Probability and Statistics; ‘ 
(5) Mathematical Physics and Applied Mathematics; 

(6) Logic and Foundations; 

(7) Philosophy, History, and Education. 


In each of these sections half-hour addresses will be delivered by experts on 
invitation of the Organizing Committee. Moreover short lectures will be given 
by members of the Congress who have applied beforehand to the Organizing 
Committee. The time allotted for each short lecture will be 15 minutes. It will 
depend on the number of these short lectures whether and how the sections will 
be divided into subsections. 

The Organizing Committee is planning several entertainments and also a 
number of interesting excursions. 

Those who wish to attend the Congress are requested to communicate their 
name (with degrees, qualifications, etc.) and full address to the secretariat 
(2d Boerhaavestraat 49, Amsterdam, Netherlands) as soon as possible. They will 
receive a more detailed communication which will be sent out in the course of 
1953. J. F. Koxsma 

Secretary of the Organizing Committee 


Heineman Foundation for Research, Educational, Charitable, and 
Scientific Purposes, Inc. 


The Board of Directors of the Heineman Foundation for Research, Educa- 
tional, Charitable, and Scientific Purposes, Inc., announce the establishment of 
a prize to be known as the ‘‘Dannie Heineman Prize’’, in the sum of $5,000.— 
to be awarded every three years to the author of an outstanding book or manu- 
script in the mathematical or physical sciences. The object of the prize is to 
encourage the writing of books on a high scientific level which have merits of 
exposition and which are likely to facilitate access to important fields of research. 

Any author desiring to submit a book or manuscipt for consideration shall 
deliver two copies to the Secretary of the Foundation. Submissions for the next 
award shall be made not later than the 31st December 1955. Detailed rules 


governing the awards may be obtained from the Secretary of the Foundation, 
50 Broadway, New York 4, N. Be UE SINE 


Vol. IV, 1953 225 


SSS Se eee 
Buchbesprechungen — Book Reviews — Notices bibliographiques 


Servomechanisms and Regulating System Design, Volume I. By 
H. CHEsTNUT and R. W. Maver (John Wiley & Sons, Inc., New York, 1951) 
#97 pp.) 350 Ties. $7.75: 


Der vorliegende Band ist vor allem ein Lehrbuch, das sich zum Ziele gesetzt 
hat, auch dem theoretisch weniger begabten Physiker und Ingenieur die Ein- 
arbeitung und rechnerische Erfassung seiner Regulierprobleme zu erméglichen. 
Dieses Ziel diirften die Verfasser erreicht haben, was hauptsachlich der durch- 
dachten Gliederung und der subtilen Darstellung des Stoffes zu verdanken ist. 
Nach einer Einfiihrung in die Problemstellung (1. Kapitel), einer kurzen Dar- 
stellung der komplexen Rechnungsweise (2. Kapitel), der Lésung von linearen 
Differentialgleichungen mit erzwungenen Vorgangen (3. Kapitel), ihrer Behand- 
lung mit der Laplace-Transformation (4. Kapitel) und dem stationaren Verhalten 
bei sinusférmiger Antriebsfunktion (5. Kapitel) beginnt die eigentliche Behand- 
lung des Stoffes auf Seite 124 mit der Diskussion der Methoden fiir die Feststel- 
lung der Stabilitat von geschlossenen Regelkreisen (Kapitel 6). Routus Krite- 
rium zur Abzaéhlung von Wurzeln mit positiven Realteilen einer charakteristi- 
schen Gleichung hohern Grades wird ohne Beweis angefiihrt. Das Nyquist- 
Stabilitatskriterium findet eine eingehende Wiirdigung und wird an Hand von 
Beispielen erlautert. Das Kapitel 7 behandelt oft verwendete Reglerteile und 
ihre Darstellung durch komplexe Ubertragungsfunktionen (Transfer Functions). 
Drei verschiedene Arten Reguliersysteme werden im 8. Kapitel postuliert, namlich 
Typ 0. Ein konstanter Wert der geregelten Grésse verlangt ein konstantes 
Fehlersignal. 

Typ 1. Ein konstantes Anderungsverhaltnis (Geschwindigkeit) der geregelten 
Grésse verlangt ein konstantes Fehlersignal. 

Typ 2. Eine konstante Beschleunigung der geregelten Grosse verlangt ein kon- 
stantes Fehlersignal. 

Nur diese drei Typen werden in der Folge betrachtet, da andere praktisch 
nicht in Frage kommen. Starkes Gewicht wird auf die allgemeine Bezeichnungs- 
weise und die Buchstabensymbole gelegt. Beides ist gemass den bereits vorhan- 
denen Empfehlungen des A.I.E.E. (American Institute of Electrical Engineering) 
konsequent im Buche durchgefiihrt worden. Die Symbole sind entsprechend der 
Funktion im Regelkreis festgelegt und kénnen je nach Aufbau ganz verschiedene 
physikalische Gréssen reprasentieren. Diese Darstellungsweise hat den Vorteil, 
dass sie in allen Gebieten der Wissenschaft Giltigkeit hat. Im Kapitel 9 und 10 
wird das Verhalten einfacher geschlossener Regelkreise mit Hilfe der Ortskurve 
der Ubertragungsfunktion in der komplexen Ebene (Nyquist-Diagramm) be- 
schrieben und verschiedene Stabilisierungsmassnahmen diskutiert. Wesentlich 
wichtiger fiir die praktische Berechnung und deshalb in den Kapiteln 11 und 
12 eingehender behandelt ist das Dampfungsverfahren, welches das Nyquist- 
Kriterium mit dem Dampfungs- und Phasenverlauf der Ubertragungsfunktion 
in Abhangigkeit von der Frequenz verkniipft. Bopes Theoreme tiber den Zusam- 
menhang von Real- und Imaginarteil von Ubertragungsfunktionen werden tiber- 
nommen und zweckentsprechend angewandt. Komplizierte Ubertragungsfunk- 
tionen lassen sich dabei sehr einfach darstellen (Dampfung als gebrochenen 
Streckenzug mit verschiedenen Neigungen von 0, 20, 40, 60, ...db/Dekade). 
Besondere Netzkarten («Nichols» Charts) erméglichen das Bestimmen des Fre- 


, 
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quenzverhaltens des geschlossenen Regelkreises. Besonders einfach gestaltet sich 
die Diskussion stabilisierender Massnahmen’ mit Zusatzvierpolen. Regelkreisen 
mit mehrfachen Riickfiihrungen und mehreren Eingangen ist das Kapitel 13 
gewidmet, wobei sich das Dampfungsvertahren zur Behandlung besonders be- 
wahrt. Endlich wird.im Kapitel 14 der Zusammenhang von Frequenzgang und 
Einschwingvorgang bei Anlegung des Einheitssprunges untersucht. Wertvolles, 
errechnetes Material in zahlreichen Charakteristiken zusammengestellt, gestattet 
anndherungsweise den Ubergang von entsprechenden charakteristischen Grossen, 
wie Uberschwingen, Zeitdauer bis innerhalb 5% des Endwertes usf., auf die dazu 
notwendige Ubertragungsdampfungsfunktion oder umgekehrt. 

Ein Literaturverzeichnis mit 114 Arbeiten, auf welche im Text Bezug genom- 
men wurde, ein Anhang von 50 Seiten mit Problemstellungen den Kapiteln zu- 
geordnet und ein Schlagwortverzeichnis vervollstandigen das Buch. 

Abgesehen von einigen unwesentlichen Druckfehlern bleibt kaum ein Wunsch 
offen, es sei denn, dass Probleme, wie Stéreinfliisse, Nichtlinearitat, Sattigungs- 
erscheinungen, unstetige Regelkreise usf., im folgenden Band mit, der gleichen 
liebevollen Darstellung ihre Behandlung finden werden. . 

Das vorliegende Buch kann jedem, der sich mit Regelproblemen und ihrer 
Berechnungsweise vertraut machen will, warm empfohlen werden. Aber auch der 
Fachmann wird gern zu dieser geschlossenen Darstellung greifen und ihr manche 
wertvolle Anregung entnehmen kénnen. H. Weber 


Ubungen zur projektiven Geometrie. Von H. HERRMANN (Verlag Birk- 
hauser, Basel 1952), 168 S., 94 Abb. und 4 Raumbilder (Anaglyphen) ; sFr. 14.-. 


Die Vertrautheit mit den Begriffsbildungen der projektiven Geometrie ist eine 
unerlassliche Voraussetzung fiir das Studium irgendeines Zweiges der Geometrie}; 
die grosse Zahl von Lehrbiichern iiber projektive Geometrie wird dadurch ver- 
standlich. Die Ubungen zur projektiven Geometrie von HERRMANN erganzen in 
erster Linie die in derselben Sammlung (Lehrbiicher und Monographien aus dem 
Gebiete der exakten Wissenschaften) erschienene Projektive Geometrie von 
W. BLaAscHKE. Daneben sind sie aber auch zum Gebrauch neben Vorlesungen 
und parallel zum Studium anderer Lehrbiicher iiber projektive Geometrie gedacht. 

Die Aufgaben schliessen an eine analytische Behandlung der projektiven 
Geometrie an. Die Sammlung enthalt dementsprechend keine Konstruktionsauf- 
gaben. Zu den 297 Aufgaben sind entweder vollstandige Lésungen angegeben, 
oder aber die Lésungen sind wenigstens so weit skizziert, dass sie der Leser ohne 
grosse Miihe selbst vervollstandigen kann. 

Durch die Auswahl der Beispiele sowie durch die vorgeschlagenen Lésungs- 


wege stellt der Verfasser den Begriff der Matrix in den Vordergrund, und zwar - 


weit mehr, als dies in der projektiven Geometrie bis jetzt tiblich war. Es sei zu- 
gegeben, dass dadurch neue Gesichtspunkte zutage treten. Uns will aber scheinen, 
dass tiberall dort, wo die Matrix nicht Abkiirzung der Schreibweise, sondern 
selbst Objekt der Geometrie ist, der Geometrie etwas Zwang angetan wird. Jeden- 
falls werden die Akzente von der projektiven Geometrie auf die lineare Algebra 
verlegt. 

Die Sammlung ist getrennt in Aufgaben aus der ebenen und Aufgaben aus der 
dreidimensionalen projektiven Geometrie. 

In der Lésung der Aufgabe 86, Abschnitt a, ist eine Unkorrektheit enthalten; 
zwei ebene projektive Felder lassen sich nicht immer so in den Raum legen, dass 
sie durch eine Zentralprojektion auseinander hervorgehen. Die Formulierung der 
Aufgabe 113 «Das Erzeugnis zweier projektiver Geradenbiischel ist eine C, oder 
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as Gerade» gefallt nicht ganz und diirfte dem Anfanger Schwierigkeiten be- 
reiten. 

Das Schlusskapitel bringt in Form von Aufgaben eine nette Einfiihrung in die 
Geometrie der Konfigurationen. Vier dusserst anschauliche Raumbilder oder 
Anaglyphen (eine Anaglyphenbrille liegt bei) von wichtigen Figuren der projekti- 
ven Geometrie beschliessen das Werk. 

Die vorliegende Aufgabensammlung bildet eine willkommene Erganzung zur 
Projektiven Geometrie von BLASCHKE; wer das Bediirfnis hat, seine Kenntnisse 
durch Aufgaben zu erharten, wird sie dankbar aufnehmen. Die Auswahl der Auf- 
gaben ist sehr geeignet, zu weitergehender Beschaftigung mit den beriihrten Ge- 
bieten anzuregen. ; 

Die Ausstattung ist wie iiblich bei dieser Birkhauser-Reihe vorztglich. 

M. Jeger 


The Classical Theory of Fields. By L. Lanpav and E. Lirsuitz (Addison- 
Wesley Press, Inc., Cambridge, Mass., 1951). 354 pp., 12 figs; $7.50. 

Das vorliegende Werk ist im wesentlichen ein Kompromiss zwischen zwei 
Richtlinien: Es soll einerseits ein einfiihrendes Lehrbuch in die Theorie des 
Elektromagnetismus und der Gravitation sein, andrerseits die allgemeinen Prin- 
zipien der klassischen Feldtheorien iiberhaupt darlegen, wie sie auch anderswo 
(Elastizitatslehre, Hydrodynamik, Ubergang zur Quantentheorie der Felder) 
gebraucht werden. Das Resultat des Kompromisses ist zwiespaltig: Einerseits 
sind die allgemeinen Prinzipien nicht in zusammenhangender Form, sondern nur 
am Beispiel der beiden speziellen Felder dargelegt; andrerseits ist der ihnen ent- 
sprechende deduktive Weg fiir eine einfiihrende Darstellung nicht ganz iiber- 
zeugend. 

_ Der Aufbau des Werkes ist im wesentlichen folgender: Vorangestellt (Kapitel 
I) ist das Relativitatsprinzip, das als Grundlage fiir die ganze Behandlung dient. 
Dann (Kapitel II) wird die relativistische Mechanik von Massenpunkten deduktiv 
aus dem durch Invarianzforderungen bestimmten Wirkungsintegral abgeleitet. 
Die Nichtexistenz starrer Kérper in der Relativitatstheorie fiihrt zum Schlusse, 
die Elementarteilchen seien Massenpunkte; infolgedessen wird im ganzen Buch 
nur mit Punktladungen gearbeitet, die Elektrodynamik kontinuierlicher La- 
dungen und Medien also nicht behandelt. 

Kapitel III bis IX entwickeln in dem so gezeichneten Rahmen die Theorie 
des elektromagnetischen Feldes in Wechselwirkung mit Ladungen (inklusive 
Optik), wieder deduktiv von einem postulierten Variationsprinzip ausgehend. 
Bemerkenswert sind hier insbesondere die Kapitel iiber Feld und Ausstrahlung 
bewegter Ladungen, in welchen viele Probleme in eleganter Weise besprochen 
werden, welche sonst kaum in Lehrbiichern angetroffen werden; so etwa die 
Wechselwirkung bewegter Teilchen bis zur zweiten Ordnung in den Geschwindig- 
keiten (nach DARwIN), die Ausstrahlung des Zweikérpersystems mit Coulomb- 
Wechselwirkung, Lichtemission durch Teilchen im homogenen Magnetfeld (Beta- 
tronstrahlung) usw. 

Kapitel X und XI behandeln die allgemeine Relativitatstheorie samt dem 
notigen mathematischen Werkzeug, dey Tensoranalysis. Als Anwendungen wer- 
den neben der Schwarzschildschen Lésung fiir das zentralsymmetrische, zeitlich 
konstante Gravitationsfeld die Friedmannschen Lésungen gegeben, welche einer 
raumlich isotropen und homogenen, zeitlich expandierenden Welt entsprechen. 

Das tibersichtlich und mit grosser Eleganz geschriebene Werk wird, trotz der 
gedusserten Vorbehalte, dank seinem nicht dem Ublichen entsprechenden Ge- 
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sichtspunkt, der Fiille behandelter Probleme und insbesondere der grossen Zahl 
schéner, mit Lésungen versehener Ubungsaufgaben jedem von grossem Nutzen 
sein, der sich in die Theorie der beiden klassischen Felder zu vertiefen wiinscht. 

M. R. Schafroth 


Crystal Growth. By H. E. Buckrry (John Wiley & Sons, Inc., New Yore 
1951). 571 pp., 169 figs.; $9.00. 

H. E. Bucktey hat sich in vielen schénen Arbeiten mit dem Problem des 
Kristallwachstums, insbesondere der Habitus- und Trachtbeeinflussung, befasst. 
Heute beschaftigen sich manche Laboratorien mit Avistallztichtung, einer Experi- 
mentalwissenschaft, die von den Mineralogen, denen die, mannigfachen morpho- 
logischen Aspekte der Mineralien eine Fiille von Problemstellungen darboten, 
inauguriert wurde. Die moderne Kristallphysik und technische Kristallkunde 
aber benotigt Kristalle bestimmter Eigenschaften; sie muss diese in geforderter 
Grosse und Reinheit synthetisieren. Bei diesen Versuchen stellen sich oft mancher- 
lei Uberraschungen ein. Wohl geben Theorien iiber den Wachstumsvorgang und 
den Einfluss der Konstitution der Lésungen auf den Wachstumsprozess wertvolle 
Anhaltspunkte, allein es liegt in der Natur der Vorgange, dass eine Beherrschung 
der Gesamtheit der Autbauprozesse noch unméglich ist. Man muss oft die Ver- 
suchsbedingungen variieren und ist dann froh, Erfahrungen anderer Forscher mit 
beniitzen zu kénnen, selbst wenn es sich nur um empirische Kenntnisse handelt, 
die der theoretischen Deutung noch Schwierigkeiten bereiten. Daher ist ein Buch 
wie das von BUCKLEY ein unbedingtes Erfordernis geworden, da es tiber viele 
solche Beobachtungen ausgezeichnet orientiert. Es wird bei der Anzeige dieses 
Werkes am besten sein, nicht auf Spezialfragen einzugehen, sondern eine Uber- 
sicht tiber den Inhalt zu geben. 

Nach einer kurzen Einfiihrung iiber Lésungen, Léslichkeit, Ubersattigung 
und die Methoden der Kristallziichtung folgen vier Abschnitte tiber die Theorie 
des Kristallwachstums, die den Zweck verfolgen, itiber verschiedene Ansatze zu 
orientieren, dieses Problem generell zu behandeln. Ideal- und Realkristall, Typen 
von Kristallisationsvorgingen, Lésungs- und Atzvorgange werden in drei wei- 
teren Abschnitten erGrtert. Fiir den von der physikalischen Seite an das Problem 
herantretenden Forscher sind die Erfahrungen besonders wichtig, die iiber die 
Habitus- und Trachtbeeinflussungen, den Einfluss von Verunreinigungen, iiber 
spezielle Wachstumserscheinungen und iiber den Einfluss der Bedingungsande- 
rungen wahrend der Kristallisation gesammelt wurden. 

Jeder Abschnitt enthalt ausgewahlte Literaturverzeichnisse, die dem Leser 
ermOglichen, die wichtigsten Originalarbeiten aufzusuchen. Die Fiille des zu be- 
waltigenden Stoffes war so gross, dass manche dem Mineralogen besonders am 
Herzen liegende Fragen, wie zum Beispiel die morphologische Charakterisierung 
einer Kristallart, zuriickgestellt werden mussten, da sich ja das Buch in erster 
Linie an den Kristallziichter wenden will. [hm aber wird es so mannigfache Anre- 
gungen und Hinweise iiber beim Kristallwachstum auftretende, manchmal er- 
winschte, oft aber unerwiinschte Erscheinungen darbieten, dass es fiir ihn 
unentbehrlich ist. P. Niggli 


Vektor- und Dyadenrechnung fiir Physiker und Techniker. Von 
E. Lour (W. de Gruyter & Co., Berlin 1950). 488 S., 34 Abb.; DM 24.-. 

Ein Hauptcharakteristikum dieses Buches besteht fraglos in seiner ausserge- 
wohnlich reichhaltigen Dotierung mit instruktiven und interessanten Anwen- 
dungen. Sie beziehen sich auf fast alle Disziplinen der Physik: Mechanik der 
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Punktsysteme und der starren Korper, Elastizitatstheorie, Hydrodynamik, Optik 
elektro-magnetische Theorie, Quantentheorie usw. Besondere Abschnitte sind 
rein mathematischen Gegenstanden gewidmet (Geometrie, Funktionentheorie). 
In diesen Anwendungen scheint mir der hervorragende Vorzug des Werkes zu 
legen. Sie nehmen nahezu die Halfte seines Volumens ein. Der angehende Phy- 
siker wird aus ihrem Studium ohne Zweifel vielgestaltigen Nutzen ziehen. Uner- 
lasslich dafiir ist allerdings eine intime Vertrautheit mit der in den beiden ersten 
Teilen entwickelten mathematischen Begriffswelt (Vektoren, Dyaden usw.), ins- 
besondere auch mit dem verwendeten symbolischen Formalismus. Im ersten Teil 
(etwa 100 Seiten) wird die Algebra dieser Gréssen dargestellt, worauf in einem 
etwas umfangreichern zweiten Teil die Analysis (Theorie der Felder) folgt. Die 
knappe Darstellung, besonders bei vielen Anwendungen — die reiche Fiille ware 
anders kaum unterzubringen —, heischt vom Studierenden intensive Mitarbeit. 
Der Verfasser benutzt konsequent die Symbolik von Gipps und JAUMANN, 
welche nach seiner Erfahrung «die klarste und verallgemeinerungsfahigste» ist. 
Der Referent bedauert lebhaft, dass es ihm nicht gelingen will, sich dieser An- 
sicht anzuschliessen, dies in erster Linie angesichts des héchst zweckmassigen 
Formalismus, welcher in Algebra und Analysis fiir mathematische Gegenstande 
der selben Natur in umfassender Weise entwickelt worden und heute allgemein 
iiblich ist. F. Baebler 


The Design and Analysis of Experiments. By O. KemprHorne (John 
Wiley & Sons, New York, 1952). 631 pp., 27 figs.; $8.50. 

Das Entwerfen von Versuchsplanen ist durch R. A. FisHEeR in den Rang einer 
wissenschaftlichen Methode erhoben worden. Diese Verfahren sind einer der 
wichtigsten Beitrage, welche die neuere mathematische Statistik zum Ausbau 
der Erfahrungswissenschaften beigesteuert hat. 

KEMPTHORNE setzt sich in diesem Werke zum Ziel, das Planen von Versuchen 
von einem méglichst allgemeinen Standpunkt aus zu behandeln und es in den 
Rahmen der grundlegenden mathematisch-statistischen Theorie zu stellen. 

Nacb einleitenden Bemerkungen tiber die Bedeutung der statistischen Metho- 
den im Bereiche der Erfahrungswissenschaften gibt KEMPTHORNE eine geschlos- 
sene Darstellung der Theorie der Mehrfachregression mit Anwendung auf die 
Streuungszerlegung, wie sie in dieser Vollstandigkeit und Strenge meines Wissens 
noch nirgends zu finden ist. Nachdem er die Wichtigkeit der zufalligen Zuteilung 
der Versuchseinheiten auf die verschiedenen Verfahren gebtihrend hervorgehoben 
hat, wendet er sich den einfachen Versuchsplanen zu — der zufalligen Anordnung 
in Blécken und dem lateinischen Quadrat. Hierauf werden die Versuche mit meh- 
reren Faktoren sehr eingehend besprochen, und zwar in allen ihren Abwandlun- 
gen. Auch die von F. YATES zuerst vorgeschlagenen Versuchsanordnungen in 
«unvollstandigen Blécken» werden einlasslich erortert. 

Wahrend Cox und Cocuran in ihrem Werk Experimental Designs!) sich eher 
an die Versuchsansteller wandten, geht KEMPTHORNE mehr auf die theoretischen 
Grundlagen ein. Dementsprechend erlautern Cox und CocurRan die Verfahren an 
zahlreichen Beispielen, wahrend sich KemMpTHorRNE mit der Angabe der den 
Planen zugrunde liegenden Formeln begnigt. 

Das ausgezeichnete Werk von KremprHorne ist vor allem fiir den beratenden 
Statistiker wertvoll; es ist aber auch allen jenen zu empfehlen, die sich der Be- 
deutung einer zweckmassigen Anlage von Versuchen in der reinen und ange- 


wandten Forschung bewusst sind. A, Linder 


1) Siehe die Besprechung in ZAMP 2, 53 (1951). 
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Hochfrequenztechnik und Weltraumfahrt. Herausgegeben im Auftrage 
der Gesellschaft fiir Weltraumforschung e. V. von Dr. Ing. R. MERTEN (Verlag 
S; Hirzel, Zurich 1951):116 5.4, 65tA Dbssh rai 


Das von R. MERTEN im Auftrage der Gesellschaft fiir Weltraumforschung 
herausgegebene Heft Hochfrequenztechnik und Weltraumfahrt enthalt acht nach 
Inhalt wie Darstellung sehr bemerkenswerte Vortrage, welche im Januar 1951 im 
Rahmen der 4. Jahreshauptversammlung der genannten Gesellschaft herausge- 
geben wurden. Nach Autoren und Themen aufgefiihrt, handelt es sich um: W. DIE- 
MINGER, Weltvaumfahrt und Ionosphadve; H. DOR1NG, Stand der Zentimeterwellen- 
technik; W. STEPP, Die Reichweite von Funkmessgeraten; G. ULBRICHT, Funknavi- 
gation mit Zentimeterwellen; F. W. GUNDLACH, Grundsdtzliches tiber Antennen fiir 
Raumfahrzeuge; KtRSCHSTEIN, Die Steuerung von Raumschiffen und thre Stabilitat; 
R. MERTEN, Funkverbindungen mit der Aussenstation; R. Moscu, Geschwindig- 
keitsmessungen nach dem Dopplerprinzip und thre Anwendung fiir Flugweitensteu- 
evungen und Bahnvermessungen. In ihrer Gesamtheit ergeben diese Beitrage eine 
wertvolle Ubersicht iiber die elektrischen Probleme der Fernsteuerung von Rake- 
ten bzw. von Raumschiffen, allerdings indem oft auf Einzelheiten der theoreti- 
schen Behandlung oder nahere Angaben von Schaltungen verzichtet wird und 
die Betrachtung vor allem auf die grundsatzlichen Fragen gerichtet ist. Viel Wis- 
senswertes entstammt friiheren Forschungsarbeiten tiber die V-2-Technik. Bei 
der weitgehenden Spezialisierung der heutigen technischen Wissenschaft ist es 
zweifellos ein verdienstliches Unternehmen, wenn — wie in der vorliegenden Ver- 
offentlichung — die Méglichkeit geboten wird, die Gesamtheit der Teilaufgaben 
in ihrer ausserordentlichen Vielgestaltigkeit kennenzulernen, welche mit moder- 
nen technischen Problemen — hier mit der Fernsteuerung von Raumfahrzeugen 
auf grésste Distanzen — verkniipft sind. Solche Probleme kénnen iiberhaupt nur 
auf der Grundlage der Spezialisierung und bestorganisierter Zusammenarbeit 
gelést werden. Das kleine, gut ausgestattete Werk sei Ingenieuren und Physikern, 
welche sich fiir die hochfrequenztechnischen Fragen der Weltraumfahrt interes- 
sieren, bestens empfohlen. F. Tank 


Advances in Electronics. Vol. III, Edited by L. Marton (Academic Press, 
New York, 1951). 357 pp., 118 figs. ; $7.50. 


Der dritte Band der von L. Marton herausgegebenen Biicherreihe Advances 
in Electyonics enthalt wiederum eine Reihe von zusammenfassenden Berichten, 
welche sich einerseits mit den Fortschritten auf dem Gebiet der Elektronenréh- 
ren, andererseits mit einzelnen Fragen prinzipieller Art aus dem Gebiete des 
Nachrichtenwesens befassen. 

Im ersten Beitrag, Field Emission Microscopy, beschreibt F. ASHwoRTH das 
Feldelektronenmikroskop und die damit gewonnenen Aufschliisse tiber Adsorp- 
tionsverhaltnisse von Fremdatomen auf Einzelkristallen von Wolfram und 
Molybdan. Das theoretisch erwartete Auflésungsvermégen ist mit etwa 2 AE 
erreicht worden, und es ist bereits méglich gewesen, den ungefahren Aufbau 
einzelner adsorbierter Molekiile abzubilden. Der weitere Beitrag Velocity Modu- 
lated Tubes von R. R. WARNECKE, M. CHoporow, P. R. GUENARD und E. L. 
GINZTON gibt eine ausgezeichnete Ubersicht iiber die verschiedenen Formen von 
geschwindigkeitsmodulierten Réhren, iiber ihre Verwendungsméglichkeiten und 
iiber die Schwierigkeiten, welche bei der Weiterentwicklung nach héheren Fre- 
quenzen und hoheren Leistungen zu tiberwinden sein werden. Wahrend die rech- 
nerische Behandlung von geschwindigkeitsmodulierten Roéhren noch relativ ein- 
fach ist, hat die theoretische Erfassung der im Magnetron herrschenden Ver- 
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haltnisse auf grosse Schwierigkeiten gefiihrt. In den beiden Beitragen von L. Brit- 
LOUIN (Electronic Theory of the Plane Magnetron) und von L. Brirtovin und 
F. Biocu (Electronic Theory of the Cylindrical Magnetvon) wird die Theorie des 
Magnetrons um ein gutes Stiick weitergefiihrt. J. E. WHITE berichtet im Beitrag 
Lube Miniaturization tiber die Entwicklung, welche zur Herabsetzung der R6oh- 
rendimensionen (Kolbendurchmesser weniger als 1 cm) gefiihrt haben. Die durch 
die Gréssenbeschrankung neu auftauchenden Anforderungen an Material und 
Konstruktion werden besprochen. Die fabrikationsmiassige Herstellung solcher 
Roéhren hat es erméglicht, bei Verwendung passender Einzelteile und neuer 
Verdrahtungsmethoden auch die entsprechenden Gerdte selber in ihrer Grésse 
stark herabzusetzen. Dies wird im Beitrag von G. SHAPIRO, Subminiaturization 
Lechniques, an Hand von zahlreichen Beispielen gezeigt. In Principles of Pulse 
Code Modulation bespricht H. F. Maver die neueste Form der heute bekannten 
Modulationsarten, welche — allerdings auf Kosten einer grésseren Bandbreite — 
eine erstaunliche Stérbefreiung erméglicht. Der Beitrag von E. A. GUILLEMIN, 
A Summary of Modern Methods of Network Synthésis, beschaftigt sich mit der 
Berechnung von Netzwerken mit vorgegebenen Eigenschaften, wahrend M. LEt- 
FER und W. F. SCHREIBER im Beitrag Communication Theory allgemeine Fragen 
prinzipieller Art besprechen, deren Behandlung in den letzten Jahren zu einer 
neuen Betrachtungsweise tiber das Wesen und die Méglichkeiten der Nachrichten- 
ubermittlung gefiihrt haben. 

Auch in dem vorliegenden Band hat es der Herausgeber verstanden, durch 
geeignete Wahl seiner Mitarbeiter und durch Koordinierung der einzelnen Bei- 
trage ein Werk zu schaffen, welches jeder, der sich tiber die Fortschritte der 
Elektronik orientieren will, mit Gewinn zur Hand nehmen wird. 

Einzig zur Aufnahme der beiden Arbeiten tiber die Theorie des Magnetrons 
méchte der Referent einen Vorbehalt machen. So wertvoll diese beiden Publika- 
tionen an sich auch sind, so stellen sie doch zwei ins einzelne gehende Original- 
arbeiten dar, welche in eine Fachzeitschrift gehoren. Die Aufnahme derartiger 
Beitrage bringt aber die Gefahr mit sich, dass die «Advances» ihrer eigentlichen 
Zweckbestimmung entfremdet werden. Und das miissten die zahlreichen Freunde; 
die sich die «Advances» in der kurzen Zeit seit Beginn ihres Erscheinens erworben 
haben, sehr bedauern. : A. A. Rusterholz 


Helicopter Analysis. By A. NixorsKy (John Wiley & Sons, New York, 
_ 1951). 340 pp.; $7.50. 

Das vorliegende Buch ist aus einer zweisemestrigen, fiir «graduate students» 
an der Princeton University gehaltenen Vorlesung hervorgegangen und ist auch 
seinem ganzen Aufbau nach vor allem als Lehrbuch gedacht. Dieser Bestimmung 
entsprechen auch die zahlreichen ausgerechneten Beispiele, die Aufgaben und 
der reiche Literaturhinweis. 

Zum leichteren Verstandnis sind die behandelten Probleme in einen ersten, 
einfacheren, und einen zweiten, mathematisch schwierigeren Teil aufgeteilt. 

Nach einer kurzen Zusammenfassung der Luftschraubentheorie behandelt der 
erste Teil den stationdren vertikalen und Vorwartsflug des Helikopters mit einer 
naheren Untersuchung der Bewegung der Rotorblatter. 

Der zweite Teil ist eingeleitet durch das Kapitel: «Mathematical Analysis of 
Problems in Dynamic Stability», welches eine umfassende Darstellung. des zur 
Lésung dynamischer Stabilitatsprobleme notwendigen mathematischen Appa- 
rates bringt. Es handelt sich hierbei um die Auflosung von Systemen linearer 
Differentialgleichungen, welche der Autor sowohl durch die tibliche Methode als 
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auch mit Hilfe der Laplace-Transformation durchfiihrt. Dieses Kapitel stellt eine 
willkommene Erganzung der vorhandenen Literatur tiber die dynamische Stabi- 
litat von Flugzeugen dar. 

Im weiteren werden im zweiten Teil die gewonnenen Resultate auf die Unter- 
suchung der dynamischen Stabilitat des Helikopters in der Nahe des Schwebe- 
fluges und im Vorwartsflug sowie auf die Berechnung der Reaktion des Heliko- 
pters auf die Steuerkrafte angewendet. 

Der Anhang enthalt Resultate einer vertieften Untersuchung der gleichen 
Probleme mit Beriicksichtigung der zweiten Harmonischen, Verwindung des 
Rotorblattes und dreieckférmiger Verteilung der induzierten Geschwindigkeit. 

Eine Zusammenstellung von verallgemeinerten Gleichungen zur Lésung ver- 
schiedener Arten von transformierten Funktionen der Steuerreaktionen, genaue 
Beschreibungen von Rechengangen und die Aufstellung von Rechenschablonen, 
die auch ungelernte Arbeitskrafte auswerten kénnen, machen die Helicopter 
Analysis auch fiir die Praxis ausserst wertvoll. K. Remenyik 


Mathematische Maschinen und Instrumente. Von F. A. WILLERS 
(Akademie-Verlag, Berlin 1951). 318 S., 258 Abb.; DM 34.-. 


Das Buch gibt eine Einfiihrung in das Gebiet der heute verwendeten Rechen- 
hilfsmittel, wobei neben den klassischen Geraten auch moderne elektrische und 
elektronische Vorrichtungen beriicksichtigt sind. Nach einem kurzen Abschnitt 
iiber Rechenschieber folgt ein Kapitel tiber Handrechenmaschinen, in welchem 
die verschiedenen Systeme sowie zahlreiche Ausfitihrungsformen beschrieben 
sind. Das folgende Kapitel itiber programmgesteuerte Rechengerate gibt eine 
Ubersicht iiber mehrere im Betrieb befindliche Rechenautomaten, eine Beschrei- 
bung der heute tiblichen Speicherverfahren sowie Hinweise auf die Rechenplan- 
fertigung. Ein Kapitel ist den Vorrichtungen zum Zeichnen und Messen von 
Kurven gewidmet, ein weiteres den Planimetern, deren viele verschiedene Aus- 
fiihrungsformen erlautert sind. Das nachste Kapitel beschreibt harmonische 
Analysatoren, das folgende Integraphen und Integratoren; schliesslich folgt eine 
Erlauterung der Differentialgleichungsmaschinen (Integrieranlagen), wobei die 
Einzelteile beschrieben und Schaltungsbeispiele gegeben sind. 

Durchweg stellt man eine tiberaus sorgfaltige Ausarbeitung des Stoffes fest; 
die vielen Photographien und Zeichnungen tragen zur Anschaulichkeit bei. Das 
Werk wird vervollstandigt durch ein Literaturverzeichnis mit gegen 900 Titeln und 
durch ein alphabetisches Sachregister. Es kann als Lehrbuch in Studium und 
Praxis wie auch als Nachschlagewerk bestens empfohlen werden. A. P. Speiser 


Review of Electronic Digital Computers (American Institute of Electrical 
Engineers, New York 1952). 114 pp., 103 Figs.; $3.50. 


Diese Schrift gibt die 19 Vortrage wieder, welche an der « Joint AIEE-IRE 
Conference» in Philadelphia im Dezember 1951 gehalten wurden und welche sich 
mit den konstruktiven und betrieblichen Gesichtspunkten programmgesteuerter 
Rechenmaschinen befassten. Insgesamt sind 11 Rechenmaschinen beschrieben, 
und zwar ausschliesslich solche, die sich in praktischem Betrieb bereits bewahrt 
haben. Die Beschreibungen selbst sind summarisch gehalten. Dafiir sind ausfiihr- 
liche Daten in bezug auf Betriebssicherheit, Fehlerquellen und deren Behebung 
gegeben ; dadurch unterscheidet sich die Schrift wohltuend von den meisten andern 
Publikationen ihrer Art und stellt fiir den Praktiker eine der niitzlichsten ein- 
schlagigen Neuerscheinungen dieses Jahres dar. A. P. Speiser 


EE 


Neuerscheinung | 


ricHTER Elektrische Maschinen 


Band II Synchronmaschinen und Einanker- 
umformer 


Berichtigte und ergainzte Neuauflage. 735 Seiten mit 
525 Figuren. Ganzleinen Fr. 46.80. 


Inhalt: I. Grundbegriffe der Wechselstromerscheinun- 
gen in Maschinen — IT. Die Synchronmaschine — II. Ein- 
ankerumformer. 


«Das Buch zeichnet sich vor allem durch eine vor keiner Schwierigkeit zu- 
riickschreckenden Genauigkeit aus, welche jetzt, wo Einzelleistung und Bau- 
stoffbeanspruchung der Maschinen den erreichbaren Grenzen zustreben, 
eine immer zunehmende Berechtigung erhalt. Es wird daher den bauenden 
Techniker bei Behandlung derartiger Grenzprobleme vor unliebsamen Uber- 
raschungen bewahren.»  . (PrmavVvEst, Elektr. u. Maschinenbau.) 


1951 erschien | 


Band | Allgemeine Berechnungselemente 
Die Gleichstrommaschinen 


Berichtigte und erganzte Neuauflage. 
XII und 630 Seiten mit 453 Fig. Ganzleinen Fr. 49.40. 


«Das Buch bedarf keiner besonderen Empfehlung, es soll sowohl dem schaf- 
fenden Ingenieur in der Praxis als auch dem Studierenden an Hochschulen 
in seiner klaren Sprache jenes Fachwissen vermitteln, auf welchem weitere 


Fortschritte gewonnen werden kénnen.» 
(A. GraBner, Osterr. Z. Elektrizitatswirtsch.) 


Band III: Die Transformatoren, und Band IV: Die I nduktionsmaschinen, 
befinden sich in Vorbereitung. (Band V erschien 1950 im Springer-Verlag.) 
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Werke der reinen ‘und aabewandtedt Mathematik: 


Mathematik 
BiEBERBACH, L.: Theorie der geometrischen Konstruktionen. . . geb. Fr. 18.70 


CaraATHEopory, C.: Funktionentheorie in zwei Banden Band I geb. Fr. 39.50 LS 


Band II geb. Fr. 27.05 << Uns 


CaratHtopory, C.: Mass und Integral und ihre Algebraisierung, erscheint 1953. 
Doetscu, G.: Handbuch der Laplace-Transformation. Band I: Die theoretischen 
Grundlagen der Laplace-Transformation .. . ot... Ygeb. Fr. 83:20 


Finster, P.: Kurven und Flachen in allgemeinen Ravuien . .. ‘geb. Fr. 15.408 


FuETER, R.: Analytische Geometrie der Ebene und des Raumes ._ geb. Fr. 24.95 


Herrmann, H.: Ubungen zur projektiven Geometrie . . ... . geb. Fr. 17.— 


Hurwitz,A.: Mathematische Werke,2 Bande ....... . geb. je Fr. 54.10 
Linper,A.: Statistische Methoden, 2. Auflage 1951. . . . . . . geb. Fr. 31.20 
Locnuer-Ernst,L.: Differential- und Integralrechnung . . . . geb. Fr. 49.90 


LocuEr-Ernst,L.: Einfiihrung in die freie Geometrie ebener Kory en geb. Fr. 12.50 — 


Ostrowskt, A.: Vorlesungen tiber Diff.- u. Integralrechng., BandI_geb. Fr. 36.40 
Band IIs geb. Fr. 69.70 


ScHLAFLI,L.: Gesammelte mathematische Abhandlungen, BandI_ geb. Fr. 56.15_ 
BandII geb. Fr.56.15 
STIEFEL, E.: Lehrbuch der darstellenden Geometrie. . . . geb. Fr. 29.65 


VoELKER,D., Dortscu, G.: Die zweidimensionale Laplace-Transf. geb. Fr. 47.85. 
ZWINGGI, E.: "Versicherungsmathematik x Wok, Sum os x) Ses Ra eebaerieam 


Physik und Aotronomie 


Gomsas, P.: Theorie und Lésungsmethoden des Mehrteilchenproblems der Wellen- — 


mechanik ... - .. geb. Fr. 33.30 
NIETHAMMER, TH.: Die genauen Methaden. der astronomisch-geographisehen Orts- 
bestimmung. . . - » « geb. Fr. 86.40 


VON DER PAHLEN,E.: Emfibeaeel in Adie Dynami von 'Stérnsystemen geb. Fr. 39.50 
WALDMEIER,M.: Einfiihrung in die Astrophysik ........ geb. Fr. 49.40 
WALDMEIER, M.: Tabellen zur heliographischen Ortsbestimmung ._ geb. Fr. 14.55 
WALDMEIER, M.: Die Sonnenkorona, BandI . . . . . ss ED. Frs2oze 
Kongressbericht der Internationalen Konferenz iiber Kernphysik und theoretische 

Physik, Basel 1949. 248 S. mit vielen Illustrationen. Broschiert Fr. 16.65 


Ingenteurwissenochaflen 


Bovea,E.: Giorgis rationales MKS-MaBsystem. Broschiert. . . . Fr. 28.10 
GOLDSTEIN, J.: Die, Messwandler . ip. aul, x chen a. oR ee De eo sn 
JAEGER,CuH.: Technische Hydraulik. . .. . . » « geb. Fr. 52.— 
KUHLMANN, K.: Theoretische Elektrotechnik. Band ll: Grundziige der Theorie der 
Wechsélstréme und des einphasigen Transformators. . . . . . geb. Fr. 27.55 
Band III: Grundziige der Theorie elektrischer Maschinen. . .  geb. Fr. 74.90 
Meissner, E., und Zircier,H.: Mechanik. Band I: Statik der starren, fliissigen und 
¢lastischen Vito) 0) Sis Subic sonia es Bede em SOR. Bindaeca ol 
Band II: Dynamik der starren Kérper nonald iy) Sinan. exqunatis ooh COD uk Cente 
Band III; Dynamik der Systeme er geb. Fr. 46.80 
MicHakEL, W.: Ortskurvengeometrie in der komplexen Zahlenebeiie geb. Fr. 18.50 
SANnGER,R.: Ballistische Stérungstheorie . . pe ae eo, 's/ie EOE OD, Ming Ola 
Srussi, R.: Vorlesungen tiber Baustatik, Band I cath gs «eee ROD. hog40:Gg 
RusTERHOLZ, A.: Elektronenoptik, BandI: Grundziige . . . . . . geb. Fr. 31.20 
RicuTer, R.: Elektrische Maschinen Band I= geb. Fr. 49.40 


Band Il geb. Fr. 46.80 
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Jetzt wieder lieferbar: 


OSTROWSKI 
Vorlesungen iiber Differential- und Integralrechnung 
Band I: FUNKTIONEN EINER VARIABLEN 
in zweitem, berichtigtem. Nachdruck 


Von Prof. Dr. A.’ Ostrowski, Professor an der Universitat Basel 
373 Seiten mit 42 Figuren. Ganzleinen Fr. 36.40, brosch. Fr. 32.25 


«Das Gesamtwerk wird ein grundlegendes Lehrbuch unserer Epoche sein, das den 
Studierenden in modernem Geist in die klassische Analysis einfihrt, ohne ihn mit 


den modernen Begriffsbildungen und Methoden zu tberfittern.» 
(W. Siiss, Archiv der M athematik. ) 


Band II: Differentialrechnung auf dem Gebiete mehrerer Variablen (195 1) - 
Band III: Integralrechnung auf dem Gebiete mehrerer Variablen (erscheint 1953). 
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